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ABSTRACT 


A description is given of the physical setting, facilities, and current program of Mauna Loa Observatory in 
Hawaii, with emphasis on its suitability as a site for field studies on a wide variety of phenomena. 


1. INTRODUCTION 


The new high-altitude observatory established in July 
1956 at an elevation of 11,150 feet on the slopes of Mauna 
loa, in the Hawaiian Islands, has already been described 
by Fox [2,3]. The purpose of the present paper is to re 
view the observatory’s first 2 years and to consider in 
somewhat greater detail the climatic environment and 
such aspects of the physical setting, facilities, and pro- 
grams most likely to interest those concerned with the 
possible suitability of such a mountain station for the 
study of specific phenomena within and outside the 
atmosphere. 

Mauna Loa (Great Mountain) is situated on the island 
of Hawaii, largest and southernmost of the Hawaiian 
To the north it faces its sister peak, 
Mauna Kea, across a 6,000-foot saddle. To the east, 
west, and south it falls gently to the sea from its 13,680- 
foot summit, comprising thereby more than half the entire 
island (fir. 2). Measured from its roots on the ocean floor 
18,000 foot below sea level, Meuna Loa is the earth’s 
greatest “iigle mountain mass. About 128 square miles 
of its surface lie above 10,000 feet. 

Geolov ically, Hawaii is the youngest of the 
Hawaiia) Islands. It is, in fact, still being formed [17}— 


a process to which Mauna Loa itself contributes by erupt- 
ing every few years. Consequently, its slopes have not 
been eroded into the deep valleys and sharp ridges of the 
older, more northerly islands, but retain the smooth and 
gradual symmetry (the average grade is 7 percent) of the 
fresh volcanic cone (fig. 3). The lower slopes are regions 
of heavier rainfall and are densely vegetated; higher eleva- 
tions are virtually barren wastes of dark lava. 

At about 19°30’ N. latitude, Mauna Loa is well within 
the geographic Tropics. A solitary peak, except for 
Mauna Kea, and more than 2,000 miles from the nearest 
continental land mass, it lies in the midst of a tropical 
ocean which, in that vicinity, has a mean annual temper- 
ature of 75° F., an annual range of 5° F., and a daily 
range of less than 2° F. 

The prevailing surface wind in this locality, south and 
west of the eastern Pacific high pressure cell, is the north- 
easterly trade, which has a frequency of 90 percent in 
summer and 80 percent for the year as a whole. Also 
characteristic of the region, and occurring with nearly 
the same frequency, is the trade inversion, whose average 
height of approximately 6,500 feet is only half that of 
Mauna Loa itself and corresponds to the mean timberline 
on its slopes. 
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Figure 1.—The Hawaiian Islands. 


2. HISTORY 


The potential value of a high-altitude geophysical 
observatory in the oceanic Tropics had long been realized. 
The uniqueness of Maune Loa in this respect—conferred 
by its height, insularity, mildness of climate at all eleva- 
tions, distance from sources of industrial pollution, yet 
accessibility from large cities, freedom (due to the trade 
inversion) from most of the water vapor and debris of the 
lower atmosphere, and by the climatic uniformity of the 
maritime tropical environment—was given early recog- 
nition by the First Pan Pacific Science Congress which, 
meeting in Honolulu in 1921, adopted a resolution calling 
for the establishment of a weather station at the summit. 
But it was not until 1950 that the laying down across the 
lava wastes of a cinder road traversable by motor vehicles 
made the upper reaches of the mountain readily accessible 
and an observing station there feasible. 

In December 1951, a small masonry hut was constructed 
at an elevation of 13,400 feet—one-fourth mile from, and 
280 feet below, the summit. This was designed to house 
90-day weight-driven recorders for atmospheric pressure, 
temperature, relative humidity, wind, sunshine, and pre- 





cipitation, and to provide temporary shelter for scientific 
parties (fig. 4). Building and equipment were admittedly 
crude, and were regarded only as a first step toward the 
realization of a manned mountain observatory. 

At the same time, rain gages and instrument shelter 
containing conventional barographs, hygrothermographs, 
and thermometers were placed along the slope at elevations 
of 5,100, 8,300, and 11,500 feet and, together with the 
summit station, were visited at intervals by employees of 
the U.S. Weather Bureau Office at Hilo, Hawaii. From 


the accumulating records, the first quantitative portrait 


of the climate of the upper mountain began to emerge 
[12]. 

In June 1956, in cooperation with the National Bureaw 
of Standards, the Weather Bureau erected at 11,150 feet 
a much larger structure intended to provide working 
facilities and living quarters for extended stays by scien- 
tific parties or for a permanent staff. This was named 
the Slope Unit of the Mauna Loa Observatory, and will 
be referred to in what follows as the Mauna Loa Observ- 
atory or simply as the Observatory. With the assignment! 
to Mauna Loa of an important role in the Inter 
national Geophysical Year (IGY) [4] a permanent staf 
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JANUARY 1959 
of three Weather Bureau employees was appointed in 
July 1957 to conduct its programs of meteorological and 
other geophysical observations. 


3. OBSERVATORY SITE 


Altitude: 11,150 feet (3,398 meters) 
Latitude: 19°32’ N. 

Longitude: 155°35’ W. 
Geomagnetic latitude: 19.9° N. 

The Observatory lies within a leveled 4.05-acre plot on 
the gently sloping north-northeast face of Mauna Loa. 
On every side stretches the vast barren sea of dark aa 
lava,' loosely heaped and crumbling, which composes the 
the upper mountain, with and there 
smaller patches of pahoehoe,' like pools of congealed 


2 
mantle of here 


molasses (fig. 5). 
4. ACCESSIBILITY 


From Hilo, the largest city (population 25,000) on the 
island of Hawaii, the Observatory is a 2-hour drive over 45 
miles of crushed lava road. Honolulu, a modern American 
city of 300,000 and the cultural, economic, and population 
center of the Hawaiian Islands, with research institutes, 
technical libraries, and a large university, is only an hour 
and a quarter distant by air, with frequent daily com- 
mercial flights. 

! The Hawaiian terms for the local lavas have entered into the vocabulary of volean- 


Aa resembles loosely strewn rubble; pahoehoe is billowy or ropy and, because 
Both are basic 


ology 
it is often relatively unbroken, much more highly reflective than aa, 
and have the same composition in the molten state. Differences arise during cooling, 
the aa forming at lower temperatures, smaller gas content, and more advanced crystal- 
See Gordon A. MacDonald, ‘‘Pahoehoe, Aa, and Block Lava,’’ American 
Journal of Science, vol. 251, Mar. 1951, pp. 169-191. 


lization 


“MAUNA LOA” 
13,680 


of Mauna Loa. 
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ISLAND OF HAWAII, T.H 


Figure 2.—Topographic map of Hawaii showing the Slope and 
Summit Units of the Mauna Loa Observatory, and the road from 


Hilo. 


The road up Mauna Loa gradually ascends the northeast 
face of the mountain, through a luxuriant tropical rain 
forest of epiphyte-covered trees, ferns, and sedge flourish- 
ing upon the lava. Along the route the average rainfall 


MAUNA KEA 





Hawaii, photographed from a model at the University of Hawaii, and showing the smooth slopes, gentle gradients, and symmetry 
Vertical exaggeration, 2.5 times. 
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Figure 4.—Summit Unit, Mauna Loa Observatory. Elevation 


13,400 feet (4,084 meters). 


increases rapidly from 145 inches annually at Hilo to 250 
inches between 2,500 and 3,000 feet (the island maximum 
exceeds 300 


inches) and then decreases steadily to an 


FIGURE 5. 
which largely covers the upper mountain. 
stands is plainly visible. 


Aerial view of the Observatory, looking south (upslope), taken in bright sunlight. 


Lighter patches are pahoehoe lava. 
The road to the summit leads off to the right. 
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estimated 15 inches at the summit. At about 6,000 feet 
the forest opens abruptly into fields of aa dotted with 
shrubs and scattered stands of trees. Soon these too grow 
sparse, and the road continues over a lava expanse broken 
only by gray patches of moss and lichens and an occasional 
shrub. 

The road is maintained as a public highway by the 
Territory of Hawaii only to near the 6,000-foot level, 
Above that point, deterioration of the unconsolidated sur- 
face by weather and use presents a continuing problem, 
and a good paved road to the Observatory, and eventually 
to the summit itself, remains Mauna Loa’s greatest present 
need. Still, except for short stretches where some difficulty 
in traction over the loose, hilly surface might be encoun- 
tered, the road to the Observatory is negotiable by passen- 
ger car, with due caution against overheating and, possi- 
bly, carburetor adjustment with altitude 
However, a 4-wheel drive is preferable and, for the addi- 
tional 8 miles to the summit, essential. For its own needs, 
the Observatory maintains two 4-wheel drive vehicles, one 
These 


make several round trips weekly and in rotation, leaving 


increasing 


an open truck, the other an enclosed carryall. 


one vehicle always at the Observatory for emergency use. | 





The very dark background is the aa |:va 
The leveled 4.05-acre plot on which the Observatory 
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5. FACILITIES 


The main building of the Observatory is a 20-by-40-foot 
concrete block structure, well insulated, and having a 
peaked roof of corrugated aluminum (fig. 6). It is parti- 
tioned internally into a comfortably furnished 16-by-20- 
foot living room-office (fig. 7), an instrument room, two 
bedrooms with accommodations for six, and a kitchen- 
dining room with propane hot water heater, cooking range, 
and refrigerator. 

Along its south side is a concrete slab, 15 by 45 feet, for 
mounting instruments Auxiliary structures 
include a wooden tower 8 by 8 feet with a platform 12.5 
feet above the ground; the generator shed; the Dobson 
and various other instrument housings; fuel and water 
tanks; and the anemometer mast. 

Electricity.—110-volt, single-phase, 60-cycle alternating 
current for household use and the scientific instruments is 
furnished by two recently installed 35 kw. diesel generators. 
These provide a more dependable source of electric power 
than that previously available, and are being operated 
alternately to permit preventive maintenance and thus 


outdoors. 


Fe 


— 
a 


& 6.—The Mauna Loa Observatory, looking north. 


iges and instrument shelter, (4) diesel fuel tank, (5) anemometer mast, (6) water tank, (7) main building, (8) concrete apron, (9 


From left to right: (1) generator shed, (2) solar radiation instruments, (3) 
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insure against power interruptions. 
approximately 10 kw., a substantial reserve remains to 
meet additional requirements. 
tors, with a total capacity of 600 watts, control the fre- 


Since present use is 
Three frequency regula- 


quency to the more critical circuits and to the recorders. 
Water is obtained by roof catchment drained into a 
1,000-gallon tank and augmented, if necessary, by haulage 
from Hilo. 
Workshop space is provided in the generator shed and in 
The Observatory 
well-chosen, 


the enclosed base of the wooden tower. 
is equipped with a small, but 
library, two desk calculators for data reduction, a micro- 
scope, the usual hand tools, a drill press, and a variety of 
electrical and electronic test instruments. Additional 
facilities for the purchase, fabrication, and repair of 
equipment are to be had in Hilo and Honolulu. 


technical 


Communications.—50-watt transceivers at Mauna Loa 


maintain an open 24-hour radio contact with the Hilo 
Weather Bureau station, and the feasibility of a direct 
radio link with tle Honolulu office, as well, is presently 
the 


being explored. Mobile units identical to those at 








rain 


Kiess- 


rliss spectrograph shelter (see reference [6]), (10) Dobson spectrophotometer housing, (11) instrument platform, (12) fission products 


lector. 


id below the level of the Observatory. 


In the background is Mauna Kea, about 25 miles distant. 


Intervening clouds are trade wind cumuli and lie over the saddle 
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FIGURE 7. 
south wall. 


Interior view of the Observatory’s living room-office, 
(1) sunshine panel; (2) wind, sunshine, rainfall re- 
(3) frequency regulators; (4) 50-watt transceiver; (5) mi- 
crobarograph; (6) from top to bottom, recorders for horizontal and 
normal incidence pyrheliometers, radiometer temperatures, and 
net exchange and total hemispherical radiometers; (7) recorder 


corders: 


for surface ozone. 


Observatory are mounted also in its vehicles. A powerful 
AM-FM-short television set 
tribute to the recreation of staff and visitors. Reception 


wave receiver aod a con- 


is excellent. 


6. STAFF 


The Observatory’s present complement includes the 
physicist-in-charge, two meteorological aides, and an elec- 
tronics technician. Except for the latter, staff members 
spend 6 successive days on Mauna Loa and the seventh in 
Hilo, where those who are married maintain family resi- 
dences. At least two persons remain on duty at the 
Observatory at any time. The principal responsibilities 
of the staff are to operate and service the scientific instru- 
ments, to take hourly weather observations, to reduce 
portions of the data, to cooperate in the projects of visiting 
scientists, and to do research. 

Although visitors do occasionally experience some of the 
well-known discomforts of mountain sickness, the staff’s 
acclimatization to the altitude has been complete and their 
general level of health excellent. 


7. CLIMATE AND WEATHER 


Unfortunately for so rich a subject, the scope of the 
present paper precludes more than a few remarks and some 
illustrative data concerning the weather and climate of the 
Observatory site. Additional information on the various 
and diverse climatic regimes of Mauna Loa, including 
its summit, is available in [12]. A fuller account is 
in preparation. 

As might be anticipated from the tropical maritime 
locale, the Observatory’s climate is comparatively mild 
for the altitude. Severe or violent weather is infrequent 
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and the rigors of life at Alpine stations virtually unknown, 
On the contrary, the brilliant skies and intense insolation, 
the moderate temperature and low humidity, experienced 
in so bizarre and remote a setting, induce in most visitors 
feelings of exhilaration and well-being. 

The Mountain cireulation.—Like other mountains, 
Mauna Loa imposes upon the surrounding atmosphere a 
local circulation often remarkably indifferent to the large- 
scale movement of air in the region. This circulation has 
a complex origin, being due in part to mechanical inter- 
ference by the mountain with the otherwise free flowfof 
air in its vicinity, in part to differential heating and cooling 
of its slopes relative to one another and to the free air, 
and in part to pressure gradients set up by temperature 
differences at similar altitudes between the atmospherie 
strata overlying the mountain and those more distant, 
The thermally induced components in this motion resemble 
an immense respiration—inward and upslope during the 
day, downslope and outward at night—which is reflected 
in characteristic diurnal variations in the weather. The 
wind shifts from night to day, and humidity, cloudiness, 
rainfall frequency, and turbidity tend to increase toward 
afternoon with the influx of air from lower elevations. 

Thus, a typical day at the Observatory may dawn bright 
and clear. Visibility is excellent. Peaks on other islands 
80 miles distant and more are distinguishable without 
difficulty. The trade inversion lies several thousand feet 
below, and trapped beneath it are the clouds and the bulk 
of the water vapor, dust, and haze. In the clear atmos- 
phere insolation is intense—frequently over 1.70 cal. 
cm.~? min.~' at true solar noon—and the temperature of 
the black lava and of the air rises rapidly. 

By early afternoon, moister air appears to be seeping 
upward along the mountain. The humidity increases and 
fractocumuli advance up the slopes. In the next hours 
the Observatory may be briefly enveloped in fog or light 
rain; but by evening the clouds have dissipated and the 
conditions which opened the day return. Nights are 
generally clear. Of course, individual days may vary 
widely from this regime. 

The role of the trade inversion in this process is not 
completely understood. Almost unquestionably it may 
impede, but does not necessarily prevent, vertical trans- 
port of air along the slopes; so that on Mauna Loa, as on 
mountains elsewhere, the local circulation dominates the 
diurnal aspects of climate. 

Through turbulent exchange, mountains tend also to 
take on the properties of the surrounding free air. This 
is more pronounced at summits (because of their greater 
exposure) than on the slopes and for seasonal than for 
diurnal factors. 

Thus the climate at the Observatory reflects the 
influences of the mountain circulation, the free air, the 
continentality of Mauna Loa’s own great bulk and surface 
area, and the tropical maritime surroundings. The data 
which follow combine the year-long record at the Obserya- 
tory with the temperature, humidity, and rainfall ‘or 
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TABLE 1.—Temperature (° F.) at or near Mauna Loa Observatory 











| | | awa 
| | | 


Nov 


| | 

| Jaa. | Feb Mar. | Apr. May June July | Aug. Sept. Oct. Dee. 
—_— — — — |————__|-— - 
Highest of record_. 62 63 | 60 65 | 68 70 | 70 | 69 67 68 64 67 
Mean maximum. -._........ .-| 48.4 4.6|/ 40.3] 53.6] 558/ 588 56.7 | 56.6 57.8 56. 5 53. 0 49.6 
Mean minimum... 29.6 28.9 31.3 | 33.7 | 34.5 | 37.1 35. 6 | 36.6 36.7 36.6 33.3 31.9 
Lowest of record _--. ae 22 20 20 | 2 24 24 26 29 29 29 23 22 
Mean daily range-..--------- 18.8 | 17.7 18. 6 19.9 | 21.3 21.7 21.1 20.0 21.1 19.9 19.7 17.7 
Years of record-.--- | 6 5 4 | 6 6 6 6 5 5 5 5 5 





several earlier years at two nearby stations, first at 11,500 
feet and later at 10,958 feet. 

Temperature.—The mean monthly maximum and mini- 
mum temperatures, the mean daily range, the highest and 
lowest observed, and the years of record are given in 
table 1. Although a detailed commentary cannot be 
undertaken here, the moderateness, even of the extremes, 
is noteworthy. The rapid rise of temperature in the 
spring, its more gradual decline in autumn, and the early 
maximum with a second peak in late summer, are all 
predominantly marine and free air in character. The 
diurnal range, on the other hand—large for the altitude 
and the oceanic Tropics—is attributable to the great 
expanse of the upper mountain and to the radiative 
properties of its dark lavas. 

Rainfall.—The chief year-round source of rainfall in the 
Hawaiian Islands is orographic uplift of the moist trade 
wind. During the cooler season, and in the absence of 
the trade, this is complemented by rainfall associated 
with frontal passages and cyclonic circulations, including 
the so-called Kona storms, and in the island interiors and 
uplands by convective showers. 

Observations made over a 4-year period at or near the 
Observatory are summarized in table 2A. 

The mean annual rainfall of approximately 25 inches 
appears to be consistent with extrapolations from island- 
wide isohyetal charts based on longer periods of record, 
but using only gages below about 6,000 feet. 

The record is much too short to justify reading a 
seasonal trend into it, but the intimation of a double 
maximum is at least not incompatible with the cyclonic 
rainfall contribution of winter and the greater frequency 
during the summer of local convective showers on the 
upper mountain. The wettest month of record was 
January 1956 with 12.89 inches. The variability of rain- 
fall is suggested by the fact that in 1958 the same month 
had only a trace. 

The indication in table 2B of a maximum frequency in 
late afternoon presumably reflects again the formation of 
showers within the ascending mountain currents, and is 
in contrast with the nocturnal maximum characteristic of 
a marine climate and found at lower elevations in regions 
exposed to the trades. 

\\ ind.—Continuous observations of wind at the Observ- 
atory site date only from July 1957, but appear to confirm 
Others made during short intervals over the preceding 
years. On the whole, speeds are moderate, and no 
ob\ ous topographic accelerations have been evident. 


498746—59-——2 


Strongest winds of the year were the sustained speeds of 
75 m.p.h. with gusts to 100 experienced in March 1958. 
This was not, however, a local condition, but a time of 
high wind generally, due to an unusual situation—the 
proximity of a tropical storm. Again in January 1959 a 
similar event occurred, this time producing average winds 
of 80 m.p.h. with gusts exceeding 105. 

Table 3 shows for January and July 1958 the average 
hourly speed and the frequencies of specified directions. 
July winds, in the mean, slightly exceeded those for 
January, but evidence of seasonality from so brief a record 
is questionable. On the other hand, the abrupt turnabout 
from southerly (downslope) flow at night to the upslope 
directions (west through north to east) during the day is 
unmistakable, and—as was mentioned earlier—is in 
marked contrast with the winds in the nearby free air, 
which exhibit little diurnal variation. The reversals in 
air movement, being related to the heating and cooling of 
the mountain slopes, follow shortly after sunrise and 
sunset (compare, in this respect, January with July); the 
lightest wind occurs at the times of the reversals. 

Humidity.—With the bulk of water vapor confined by 
the trade inversion to levels below the Observatory for 
most of the day, the average absolute humidity is low. 
The relative humidity displays two noteworthy charac- 
teristics: first, a diurnal course quite unlike that com- 
monly observed, in that humidity on many days tends 
to vary directly, rather than inversely, with temperature 
(fig. 8); and second, a high frequency of exceptionally 
low humidities—often below 5 percent for extended 
periods, and occasionally below 1 percent for most of a 
day. The first proceeds directly from the invasion of 
the upper slopes by lower level air during the afternoon. 
The second is more obscure, being at times almost certainly 


TaBLE 2.—Precipitation at Mauna Loa Observatory 


—oOo—[{—{{{_—_—yx<=_<__>—<_<_—_<_<_<_<_$—$—>>|>|>—_—=—_—_{——=—_={_—>[>=>[=>=_=_=_____——"— = = = 


B. Hourly Rainfall, Percentage Fre- 
quency (trace or more) of All Observa- 
tions 


A. Mean Monthly Rainfall 


Month Inches Hour ending| % || Hourending, % 
(LST) (LST) 
January 3.9 01 20 13; 49 
February 2.1 02 2.6 14 6.6 
March 4.0 6; 20 15 6.1 
April 0.7 04 20 16 6.4 
May.... 12) 05 23 17 7.5 
June... 0.3 | 06) 26 18 6.4 
July 20) 07 | 26 19 4.0 
August _- 29) os} 29 20 4.0 
September 1.9 wo; 23 21 2.6 
October ___- 1,2 | 10 | 3.2 22 17 
November - 1. 5 |} 11 3.5 23 L7 
December - 3.4 12 4.0 24 23 


! 
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TaBLe 3.—Mean hourly winds at Mauna Loa Observatory 
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January 1958, Monthly Mean 9.1 m.p.h. 


Hour ending (LST) Ol 02 03 04 O05 06 07 os 09 10 
Percentage W 

through E 19 13 16 13 16 16 16 16 19 71 
Percentage SE 

through SW 81 87 4 87 M4 M S4 MM 77 26 
Mean wind speed 

(m.p.h 9.4 9.4 9.6 9.8 10.7) 10.3 9.9 9.8 7.0 5.5 


July 1958, Monthly Mean 9.5 m.p.h. 


Hour ending (LST) 01 02 03 04 05 06 07 08 09 10 
Percentage W 

through E 0 0 0 0 0 0 0 39 81 S4 
Percentage SE 

through SW 97 97 97 100 O4 100 100 52 19 16 
Mean wind speed 

(m.p.h.) 11.5 | 10.6) 10.8) 10.3 | 10.5) 10.1 9.4 7.1 7.3 8.3 


“Indicates sunrise and sunset. 
Underscore indicates minimum wind speed. 
Double underscore indicates maximum wind speed. 


due to nocturnal subsidence over the radiatively cooled 
mountain, but at other times—as when low humidity 
persists through the day—to the breaking down of the 
ascending mountain current (perhaps when insolation is 
weak and the inversion strong), or to marked subsidence 
within the general air flow around the southern flank of 
the Pacific High. 

Measurements of the total precipitable water above 
Mauna Loa [16] suggest that the afternoon rise in surface 
humidity involves only a relatively thin skin of moister air. 


Clouds and sunshine.—Visual observations, from which 
the diurnal and seasonal variations in cloud types and 
frequencies could be obtained, are available only for a 
portion of the daylight hours and the single year that the 
Observatory has been manned. Other useful indices of 
sky cover are the proportion of possible sunshine received 
in each daylight hour and the frequency of clear, partially 
obscured, and overcast skies. Observations of the latter 
kind are summarized in table 4 for August 1957 through 
February 1958. 

Since even the thinnest and most transparent clouds 
were included in this tabulation, it considerably under- 
states the sunniness of the site and the frequency with 
which the solar dise can be seen; and it would appear 
from these and other data that about 90 percent of the 
nighttime and well over 50 percent of the daytime hours 
are nearly or entirely free from cloud above the level of 
the Observatory. 

The afternoon maximum in sky cover is due again 
primarily to the formation of fractocumuli in air rising 
along the mountain slopes. Perhaps because of its 


TaBLe 4.—Mean hourly sky cover (percentage frequency) at Mauna 
Loa Observatory 

Hour of observation (LST) 08 | 09 | 10 | 11 | 12) 13 | 14) 15) 16 17 

Clear skies 59 | 56 | 52| 47/ 38 | 29| 21/18 | 19| 24 

Half-covered sky or less 7 71 | 70 | 66 | 59 | 49 | 37 | 37 | 38 46 

Overcast 16 | 14 | 17 | 23 | 20 | 28; 31 | 35 | 31 | 25 
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smooth terrain and gentle slope, Mauna Loa has fewer 
orogenic cirrus than do some of the lower but more abrupt 
mountain barriers of the other islands. 

Snow.—Although it has on the average several snowfalls 
each winter, Mauna Loa does not even then bear a perma- 
nent snow mantle. Snow, usually light, also occurs at the 
Observatory from time to time, but the mean snow line is 
at about 12,000 feet. With intense insolation and low 
humidity, evaporation is extremely rapid; vet snow has 
been found at the summit even in midsummer, tucked 
away beneath the lava. 

Thunder, lightning, hail. 
these are infrequent and relatively mild at the Observa- 
tory. From July 1957 to August 1958 only 5 days had 
thunderstorms, and on all these occasions hail approx- 
imately one-sixteenth inch in diameter was also observed. 

“Seeing” and turbidity.—Visiting astronomers [6] have 
found the “seeing” highly satisfactory, but systematic 
observations of it, day and night, have not yet been 
attempted. However, computations of turbidity from 
solar radiation measurements (q.v.) will soon begin. 
Scattering in the circumsolar region appears to be minimal. 


8. SOME SCIENTIFIC USES 
IN THE PAST 

In the 2 years since the construction of the Observatory, 
a number of scientific endeavors have taken advantage o/ 
its elevation, the low turbidity and water vapor content 
of the overlying atmosphere, or some other special attr 
bute of the site. These studies are briefly listed below, 
with references to published reports wherever possible. 

The forms taken by snow crystals growing under natural 
conditions in an aerosol-free environment, December 1956 
January 1957, Geophysical Research Directorate [9 
With the buildings at 11,150 feet providing a base of 
operations, microphotographs were taken of freshly fallen 
snow flakes at the summit of Mauna Loa. 

Spectrographic observations of water vapor in the aimos 
phere of Mars during that planet’s close approach i 
July 1956 and of Jupiter in May 1957, National Geo 


Like other severe weather, 
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A representative hygrothermograph chart, Mauna Loa Observatory. The tendency for humidity to increase during the after- 


noon probably reflects the influx of air from lower levels. 


graphic Society and National Bureau of Standards [6]. A 
high-dispersion laboratory-type instrument was_ used. 
The prime requisite was a minimum of water vapor inthe 
terrestrial atmosphere. 

Spectroradiometry of the sun between 0.3 and 2.5 microns, 
June 1956 and May 1957, National Bureau of Standards 
16]. A lead sulfide cell and a photoemission tube were 
used in these observations to cover the spectral range. 
Strong absorption by water vapor in this band required 
that the superambient atmosphere be as free from it as 
possible. 


Atmospheric transmission in the infrared, July to Sep- 
tember 1957, Naval Research Laboratory [19]. Light 
from two 5-foot carbon-are searchlights mounted at about 
10,000 feet on the facing (south) slope of Mauna Kea was 
analyzed by an infrared spectrometer at the Observatory, 
7 miles distant. This project utilized the long optical 
path through the clear air above the inversion offered by 
the facing sister peaks (fig. 2). 

Lunar occultation program, October 1957, Army Map 


Servi Precision timing of stellar occultations by a 
photo-multiplier operating off a 12-inch reflecting telescope 


Was used to connect the Hawaiian Islands with the North 
American geodetic datum. 


PRESENT PROGRAMS 


Present observations at Mauna Loa are of three kinds: 
local weather; continuation of the International Geo- 
physical Year (IGY) programs in atmospheric composition, 
solar radiation, and particulate matter; and others for 
special purposes. 
and synoptic analyses are made also at Hilo and Hono- 
lulu.? 

Routine 
the Observatory, local weather is recorded hourly between 
0600 and 1700 Ls? and autographic tracings are obtained 
of the atmospheric pressure, temperature, relative hu- 
midity, wind speed and direction, precipitation, and 
duration of sunshine. Time-lapse cloud photography is 
used routinely to preserve the visual aspect of the weather. 

These observations serve not only the usual meteorolog- 
ical and climatological purposes, but by describing the 
local state of the atmosphere they provide the correlative 
information needed to validate and interpret data from 
the IGY and other programs. An example of this is the 
present attempt by the authors to account for the diurnal 


Pertinent meteorological observations 


meteorological observations and analyses.—At 


2 Programs in geomagnetism, solar flares, and cosmic rays begun during the IGY are in 
progress elsewhere in the Hawaiian Islands, including the 10,000-foot summit of 
Haleakala, Maui, by the University of Hawaii and other organizations. 
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recorder. Arrow 


surface 
points to glass tubing from air intake mounted 3 feet above the 
Observatory roof. 


FIGURE 9. tegener automatic ozone 


variation in surface ozone at Mauna Loa in terms of the 
mountain circulation. 

Further correlative information is provided by the 
Weather Bureau’s aerological station at Hilo where upper- 
air soundings of temperature, humidity, and wind, often 
to altitudes of 30 km., are made twice daily (at 0200 and 
1400 ist, four times daily on IGY World Days), and by 
the extensive analysis and forecasting program of the 
Honolulu Weather Bureau Airport Station, 210 miles to 
the northwest, which prepares synoptic charts of the 
broadscale circulation of the atmosphere at sea level 
(4 times daily) and at 700, 500, and 300 mb. (twice daily) 
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Figure 10.—Diurnal variation of surface ozone, November 19657, 
The afternoon minimum may be due to the influx of lower-level 
air whose ozone content has been depleted by contact with vegeta- 
tion or aerosol-bearing water droplets. 


over an area which extends from the western United States, 
through the northern and equatorial Pacific Ocean, to 
eastern Asia. 


The IGY observations..-The meteorological programs § | 
initiated at Mauna Loa during the IGY are being con- | 
ducted in accordance with the appropriate instruction & | 
manuals, and these should be referred to for further pro- & ‘ 
cedural details concerning them [15]. Instruments, proce- | 
dures, and a few preliminary results will be touched on & ! 
briefly in what follows and, wherever possible, references & ' 
will be cited for the less familiar observations. 

{ 


Surface ozone.—Continuous recordings of surface ozone 
have been made since August 1957 by means of automatic 
chemical sampling equipment recently developed by 
Regener [14]. The procedure, which is absolute and 
highly specific for ozone, utilizes the differential oxidation 
of potassium iodide to iodine within two air streams, one 
of them heated to 300° C. to dissociate its ozone. Sodium 
thiosulphate is used to prevent volatilization of the iodine. 
The resolving time of the instrument is approximately 
one minute. 

A preliminary study of the first 14 months of record 
indicates that the variations in surface ozone at Mauna 
Loa are related to both local and large-scale features of 
the atmospheric circulation. A marked seasonality 3 
evident, with a spring maximum and autumn minimum 
like that found elsewhere in surface and total ozone. 
Monthly means ranged from 33 micrograms per cubic 
meter in October 1957 to 68 micrograms per cubic meter 
in April 1958, the rise being especially steep from late 
winter to early spring. The lowest mean hourly value 
observed was about 20 and the highest near 100 micro 
grams per cubic meter; the lowest and highest 1-minute 
values were approximately 7 and 135 micrograms per 
cubic meter, respectively. 

Superimposed upon these seasonal trends, and often 
of comparable amplitude, is a pronounced diurnal varia tion 
with early morning maxima (a secondary peak is common 
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shortiv after sunrise) and minima during midafternoon 
(fig. \0). This js attributed to the local circulation of 
up- 2 d down-slope winds generated by the radiational 
heating and cooling of the mountain itself, the afternoon 


minimum being more specifically ascribed to the influx 
from below of air whose ozone content has been depleted 
during its ascent by contact with vegetation, water 
droplets, or aerosols in the moist subinversion layer 
(below about 6,500 feet). 

Striking variations in the amount of surface ozone 
have also been observed to accompany changes in the 
broad circulation patterns in the Pacific area. Thus, 
during February 1958 ozone increased abruptly with the 
transition from the prevailing trade wind regime to a 
strong westerly flow around deep cyclones centered some 
distance north of the Hawaiian Islands, and decreased 
again with the return of the trades. A report on these 
observations is in preparation. 

Total ozone-—The Observatory’s Dobson 
photometer (No. 63), by which the total depth of 
atmospheric ozone is determined from the differential 
absorption by ozone of two spectral lines in the ultraviolet 
portion of the solar beam, was put into operation on 
November 27, 1957. Its elevation (it is the highest 
permanent installation of its kind), the clarity of the 
overlying atmosphere, and the relatively small ozone 
variations typical of the latitude are expected to make 
its observations valuable not only for their contribution 
to the synoptic ozone program, but also for calibration 
principally the determination of the extra- 
terrestrial constants. The Mauna Loa spectrophotometer 
is to be used as a standard for the United States Dobson 


Spectro- 


purposes 


network. 

Lunar observations from 3 days before to 3 days after 
full moon are also being made, and these will be compared 
with those in the Antarctic where changes in the moon’s 
angular elevation are too small for close calibration of the 
instrument. 

On the basis of the first 10 months of record, day-to-day 
variations in total ozone at Mauna Loa appear to be 
small, but the amounts themselves and the amplitude of 
the seasonal variation, which ranged from approximately 
250 cm. in February to .288 cm. in May 1958 (the 
Vigroux absorption coefficients were used in these com- 
putations) are high for the latitude [13]. However, these 
values and the late occurrence of the seasonal peak 
relative to that in middle and high latitudes must await 
the confirmation of additional observations. 

Carbon dioride.—As part of the worldwide effort begun 
during the IGY to obtain baseline measurements suitable 
for defining the present CO, content of the atmosphere 
and thus detect secular trends, this component is being 
continuously recorded at Mauna Loa Observatory by 
meaiis of a Model 70 Infrared Gas Analyzer (fig. 11), an 
instrument specifically designed for the monitoring of 
flowing gas samples [1]. 

The analyzer operates by measuring the energy loss 
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Figure 11.—Infrared analyzer for CO». (1) freezer to remove water 
vapor which would otherwise interfere with the analysis, (2) an- 


alyzer box, (3) recorder. 


of an infrared beam traversing a gas sample. Radiation 
from an infrared source, after being mechanically chopped 
at 20 c.p.s., passes through the test-gas-filled cell, where 
its absorption produces a cyclic pressure pulsation. This 
is transmitted to the tantalum diaphragm of a condenser 
microphone and thence converted to a voltage which is 
recorded. A reference gas of known CO, content is used 
as a standard. Air intakes are mounted on four 20-foot 
towers 450 to 500 feet distant and in mutually per- 
pendicular directions from the main Observatory buildings. 
Two upwind flows are sampled concurrently, any difference 
between them interpreted as due to local 
contamination. 

Data secured at Mauna Loa since sampling began in 
March 1958 are sufficiently similar to those obtained else- 
where to support previous indications that in the absence 
of local pollution the carbon dioxide concentration in the 
lower atmosphere is substantially uniform over the earth. 


being 
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From 
left to right, Beckman and Whitley (Gier and Dunkle) total 
hemispherical and net exchange radiometers, and Eppley normal 
Notice the 


Figure 12.—‘Solar radiation instruments, facing northeast. 


incidence and horizontal incidence pyrheliometers. 
loose, rough texture of the ground. 


The Observatory data also suggest a diurnal variation, 
perhaps related again to the mountain circulation, and a 
seasonal trend, with mean monthly values decreasing 
gradually from approximately 314 p.p.m. (parts per mil- 
lion) in May to 312 p.p.m. in August. 

Solar The radiation measurements 
made at Mauna Loa since November 1957 are part of an 
international program designed to increase our knowledge 
of the radiation fluxes to and from the earth’s surface, and 
hence of the earth’s heat balance and budget. The in- 
struments being employed for this purpose at Mauna Loa 
are shown in figure 12. They include the Eppley (10- 
junction) horizontal incidence pyrheliometer (pyranom- 
eter) for the total incoming (short wave) radiation from 
sun and sky; the Eppley normal incidence pyrheliom- 
eter, for direct solar radiation (solar intensity), equa- 
torially mounted and driven by a synchronous motor to 
follow the sun; the Beckman and Whitley (Gier and 
Dunkle) total hemispherical radiometer for the total (long 
wave, terrestrial) radiation from earth and atmosphere; 
and the Beckman and Whitley net exchange radiometer 
for the net transfer of (long wave) radiation between earth 
and atmosphere ({15 Part VI], [5]). All are continuously 


radiation. solar 


recording. 

Filters (Nos. OG1, RG2, and RGS8) for subdividing the 
spectral regions are to be used with the normal incidence 
pyrheliometer, thus—among other things—permitting the 
turbidity coefficients to be computed. 

Fission products.—The use of the radioactive products 
of nuclear tests as tracer materials for studying the atmos- 
pheric circulation has been described by Machta [8]. As 
one of the network of fallout stations, Mauna Loa has 
operated, since August 1957, a small (40 e.f.m.) air filtra- 
tion system in which a week’s run constitutes a single 
Filters are sent for analysis to the Naval Re- 
search Laboratory, Washington, D.C. 

The Bigg-Warner expansion chamber 


sample. 


kre ng nucle i. 
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Figure 13.—Bigg-Warner expansion chamber for measuring the 


concentration of freezing nuclei. 


developed by the Australian Commonwealth Scientific and 
Industrial Research Organization for measuring the con- 


centration of freezing nuclei [18] is shown in figure 13. 


Freezing nuclei in the refrigerated air sample trigger the 
formation of ice crystals in the fog produced on expansion. 
These crystals fall into a dish of slightly supercooled sugar 
solution at the bottom of the 10-liter tank, where they en- of 
large rapidly and may be readily counted by inspection. 
By varying the initial overpressure to which the air sample 
is subjected, a spectrum of the freezing nuclei population 


th 


as a function of temperature may be obtained. 

Daily counts made at the Observatory since December 
1957 show preliminary evidence of occasional peaks oe- 
curring within a few days of those observed elsewhere. A 
current hypothesis in the field relates anomalous freezing 
nuclei activity to meteor showers a month earlier and to 
singularities in rainfall, cirrus cloudiness, and other gee 
physical phenomena. The Mauna Loa data are included 
in a recent summary by Kline and Brier [71 of these ap- 
parent correlations. 
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Micvometeorites.—The scavenging of cosmic dust from 
large volumes of air pumped through fibreglass filters 
began at Mauna Loa in August 1957 at the request of 
the Swedish oceanographer, Hans Pettersson, who saw 
in the deposition rate of micrometeoritic materials on the 
earth’s surface a possible means of dating the oceanic 
sediments [10]. After a weekly 24-hour run, the filters 
are analyzed for iron, nickel, and cobalt in the IT Chem- 
isches Institut der Universitat, Vienna. Nickel, believed 
far in excess of any possible contamination from terrestrial 
sources, has been found in most of the fiiters from Mauna 


Loa [11]. 
FOR THE FUTURE 

The year since the assignment of a permanent staff to 
the Observatory has been occupied largely with the in- 
stallation, calibration, and maintenance of the scientific 
equipment and, more recently, with the reduction of data 
for forwarding to IGY collection centers. Because of 
their value virtually all the IGY programs at Mauna Loa 
have been continued beyond the official termination of the 
International Geophysical Year in December 1958. 

As the instrumental and observational aspects of the 
work settle into routine, a limited amount of time is also 
becoming available for a preliminary study and interpre- 
tation of the record and, on a modest scale, for additional 
investigations which reflect the personal interests of the 
staff. 

Mountain circulation and climate.—The motions and 
properties of the lower 1,000 feet of Mauna Loa’s atmos- 
pheric envelope, and the structure and role in the mountain 
circulation of the trade inversion contiguous to the slopes, 
are to be explored—at first by wiresonde profiles and 
soundings of temperature, humidity, and wind, and later 
by more elaborate means. In addition to their intrinsic 
interest and their pertinence to an evaluation of the local 
factors in Mauna Loa’s weather and climate, it is hoped 
that these observations may permit the exchange coeffi- 
cients within the mountain’s envelope and between it 
and the surrounding air to be computed. 

Atmospheric interrelationships.—The ready accessibil- 
ity at the Observatory of the weather records on the one 
hand (local, upper air, synoptic charts) and the IGY 
observations on the other (ozone, carbon dioxide, solar 
radiation) makes it a convenient place in which to study 
the interrelationships between them; that is, the meteoro- 
logical correlates of diurnal, and synoptic 
Variations in ozone, carbon dioxide, and radiation, and 
conversely, the possibility of anticipating synoptic develop- 


seasonal, 


ments —for example, the formation and motion of Kona 
storms —through precursory variations in, say, total ozone. 

thers. Among other projects now in prospect for 
Meuna Loa are observations of the vertical gradients of 
ozone (by ozonesonde) and carbon dioxide (by flask 


sampling and aireraft traverses) along the mountain 


slopes and in the nearby free air. More accurate deter- 
nations of surface humidity and of the total precipitable 
Wate: in the overlying atmosphere await the receipt of 
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stable infrared hygrometers. The local and synoptic 
concomitants of the exceptionally low humidities com- 
mented upon earlier are also to be looked into. The 
freedom of the site from the pollution almost invariably 
present in continental air masses makes it well suited to 
baseline observations of atmospheric electricity, and 
measurements of conductivity, potential gradient, and 
ionization rate are being planned. 

These and other proposed studies must, of course, be 
contingent on the availability of suitable instruments and 
on prior obligations to the IGY and _ other present 
programs. 


9. CONCLUSION 


The elevation and accessibility of the Mauna Loa Ob- 
servatory, the comfort in which high-altitude work can 
be done throughout the year, its insularity, and its isola- 
tion during the night and much of the day from a large 
part of the tropospheric water vapor and turbidity, are 
advantageous for many observations within and outside 
the atmosphere. It is especially suited for those which 
require exceptional transparency of the air, and has al- 
ready been utilized in a number of such projects. 

The usefulness of the Observatory, particularly for the 
determination of standard or baseline values of, and secular 
trends in, solar radiation, the atmospheric constituents, 
and for other observations is further enhanced by its geo- 
graphical isolation. Scarcely the slightest prospect exists 
that within the foreseeable future the site can deteriorate 
significantly through the encroachment of cities or industry 
upon the present highly stable population and specialized 
agricultural economy. 

A good paved road and, less immediately, commercial 
electric power remain the Observatory’s most pressing 
needs. However, increasing recognition that the present 
site at 11,150 feet and the summit 2,500 feet above are 
potentially among the world’s most valuable terrestrial 
platforms, leaves little doubt that these improvements will 
be made. 

Present research at Mauna Loa centers about the IGY 
observational programs and on studies of the mountain 
climate and circulation, but additional work in several 
fields is planned. 

Within its limitations of time, space, and personnel, the 
staff is prepared to offer housing and its other facilities 
to visiting scientists and to assist in the instrument- 
reading-and-servicing aspects of projects established there. 
Organizations and individuals desiring further informa- 
tion about the Observatory or the use of its facilities are 
invited to address their inquiries to the Director of Me- 
teorological Research, U.S. Weather Bureau, Washington 
25, D.C., or to the authors. 
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A NOMOGRAM FOR THE DETERMINATION OF 
SOLAR ALTITUDE AND AZIMUTH’ 
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[Manuscript received August 22, 1958; revised December 3, 1958] 


ABSTRACT 


A combination nomogram is presented for the calculation of solar altitude and azimuth for any latitude, solar 
declination, and time of day. 


1. INTRODUCTION 


Meteorologists, and other scientists or engineers working 
in allied fields, frequently wish to determine the solar 
altitude and/or azimuth for various purposes. The 
present paper describes the construction of a combination 
nomogram for the calculation of both these quantities for 
any latitude, solar declination, and time of day. 


2. DERIVATION AND CONSTRUCTION OF THE SOLAR 
ALTITUDE NOMOGRAM 


The altitude of the sun is given by 


sin a=sin @ sin 6+cos ¢ cos 6 cos h (1) 
where: 
a=solar altitude (angular elevation above the 
horizon) 
¢=geographic latitude of the observer 
6=declination of the sun 
h=hour angle of the sun 
Equation (1) may be compared with the relationship 
among the coordinates of three points P;(.X,, Y1), P3(X3, Ys) 


and P;(0,¥Y;), lying in a straight line, as illustrated in 


figure 1. 
y X 
r.{ ¥! ‘_— 1 }—riZ+¥i=0 (2) 


YY, =K, sin a+C\, (: 


If we set 


~~ 
— 


Y;=K; cos h+C;, (4) 
and substitute these in (1) for sin a and cos h, we get 


Y, —p cos @ COS s\— { ‘+k, sin ¢ sin 6 


K 


— Os Fe €08 6 c08 6 }+ri=0. (: 


~' 
Sa 


If we now postulate that the left hand sides of (2) and (5) 


blished by permission of the Director, Meteorological Branch, Department of 


Transport, Canada. 


498746—59——3 


are to be identical, then 


XxX, ky 
——]=~—=> cos ¢ Cos 6 6) 
X3 K; { 
la Xi ( : K 
—'—(0,+ K;, sin ¢ sin 6——2,— cos ¢ cos 5. (7) 
*2; Ks 
In summary, the following equations define the 


nomogram 
Y, = kK, sin a+( ‘ 


Y,=K,; cos h+C; 
Xx K 


Y. — 3 cos @ cos 6 (8) 
#43 5 
eet ae ee 
Y;=( st (C\—C;+ A, sin ¢ sin 4). 
“ly 


Ky, Ci, Ks, Cs, X, are arbitrary factors, each having the 
dimension of length, which determine the size of the 
nomogram and the relative positions of the various scales. 


P, (0,y,) 
Ps (X5,Ys) 


Pi (xiy,) 











~~ X 





Figure 1.—Relationship among coordinates of three points lying 
in a straight line, similar to equation (2). 
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SOLAR DECLINATION 


Figure 2.—Nomogram for calculation of solar altitude and azimuth for any latitude, solar declination, and time of day. 


Once values are assigned to these arbitrary factors, equa- 
tions (8) may be used to calculate the positions of P, 
(sin a scale), P; (cos h scale), and P; (family of lines of 
constant ¢ and 6, hereinafter referred to as the ¢—6 grid). 
[It should be noted that the equations determining X; and 
Y; are symmetric in ¢ and 6. It follows that the line in 
the ¢—6 grid for ¢=N° is also the line for 5=N°. 

For the purposes of the original drafting, the following 
values were assigned to the arbitrary factors: K,=—20 
inches, C,=10 inches, K;=10 inches, C;=0, X;=20 inches. 
The resulting nomogram is shown in Figure 2.? 


3. USE OF THE SOLAR ALTITUDE NOMOGRAM 


The True Solar Time (also referred to as Local Apparent 
Time or Sun-Dial Time) determines a point P; on the cos h 
scale, which is graduated directly in terms of time. The 
solar declination and geographic latitude determine a 
point P; on the ¢—6 grid. If P; and P; are joined by a 
straight line, which is extended to intersect the sin a scale 


2 A working size of the nomogram, about 9 in. x 9 in., will be made available in sm ill 
quantities at no charge to interested readers. Address inquiries to M eteorological Service 
of Canada, 315 Bloor St. West, Toronto 5, Ontario. 
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at /;, the solar altitude (actually sin a) may be read off 
directly at P;. 

For use at a single station, a nomogram may be con- 
structed with only the latitude line for that station, cali- 
brated re solar declination or directly re time of year. 
However, for single station work, the graph described by 
Schutte [2] is easier to construct and more convenient to 
use. 

The solar altitude (sin a) scale may, of course, be la- 
belled directly in terms of any quantity which is a single- 
yalued function of solar altitude; e.g., optical air mass 
(average refraction assumed), CSCa at 22 km. (for ozone 
calculations), ete. 


4. DISCUSSION OF ERRORS IN THE SOLAR 
ALTITUDE NOMOGRAM 


The errors involved in the use of nomograms may be 
classified as follows: 


(a) Errors in construction (drafting, printing, warping of 
paper). 

(b) Errors in alignment (P3, Ps). 

(c) Errors in reading the answer (P;). 

(d) Curvature of straight-edge used. 


The fourth type of error may be reduced by performing the 
calculation a second time, reversing the straight-edge end 
for end in doing so, and averaging the two answers. 

The errors of types (b) and (c) result from the difficulty 
in precisely positioning the straight-edge directly over the 
points P; and P; and in reading the position of point P;. 
In general these errors will depend on the angle at which 
the straight line P,P; crosses the various scales. If we 
designate the standard error of positioning the straight 
line as A units of length, perpendicular to the line itself, 
at P; and P; and if 8 is the angle which P;P; makes with 
the X-axis, then it may be shown that the standard error 
of P; due to alignment error is, in units of sin a, 








Ase 8 (2 Ki. ) 
(1-25 RK. cos ¢ cos 6 (9) 


If we consider a further error in the final answer, due to the 


difficulty of reading off the position of P;, equal to i. a B 
1 





then the final error is 


Aeo ki +(1 —2 ie cos ¢ COS i). (10) 

For the solar altitude nomogram in a convenient work- 
ing size (about 10 in. x 10 in.), the maximum error 
occurs when ¢=6=hA=0 and is 0.724. The smallest error 
occurs when ¢=90°, 8B=0 and is 0.144. 

To provide a practical test of the accuracy obtained in 
using the nomogram, a set of 25 test computations was de- 
signed to sample all areas of the ¢—6 grid and the time 
anc! sin @ axes. These were arranged in random order 





an! given to four persons unfamiliar with the nomogram. 
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After a few minutes instruction, each person was asked to 
go through all computations once, then to repeat the entire 
set with the straight-edge reversed. The two sets of 
answers were then compared and any deviations exceed- 
ing .007 were reconciled by repeat computations. Finally, 
the reconciled answers were averaged. The standard 
error o, of the individual answer was as follows: 
Observer A (university graduate)o=0.013 (sin @ units) 


B( « “0.010 
C (high school graduate) 0.013 
a. “  ) 0.012 


If the errors in the determinations are classified accord- 
ing to geographic latitude, the following pattern of stand- 
ard error emerges 


0<¢<.20° o«=0.027 (ef. 0.72 A) 
20<o<40 o=0.012 
40<¢@ o=0.005 (ef. 0.14 A) 


Thus, neglecting errors in construction, one might infer a 
value of A=0.04 in. 

To give some idea of the accuracy obtainable by an 
experienced observer, one of the authors performed the 
same set of test computations. The standard error of 
these computations was 0.003 and a separation of errors 
according to latitude suggested a value of A=0.01 in. 
Errors of this magnitude are quite acceptable for many 
meteorological and engineering uses. 


5. SOLAR AZIMUTH NOMOGRAM 


The solar azimuth angle, a (measured from north), is 
given by 
Os 
sin a=sin h°S* (11) 
COS a 
For solution of this equation, it is convenient to use a 
proportional chart of the type described by Levens [1]. 
Such a chart is illustrated in skeleton form in figure 3. 
From pairs of similar triangles (AP;P, and BP,P,, APsP; 
and BP,Ps), it is readily seen that 


BP, BP, 


— = ° 9 
AP, AP, ae} 


Interpreting these distances in terms of the variables 
(a, h, 5, a) as illustrated, equation (12) is directly equiva- 
lent to (11). The proportional chart of figure 3 has been 
incorporated in the solar altitude nomogram as described 
below. 
The existing solar altitude scale may be used for cos a, 
noting that 
1—cos a=1—sin (4r/2—a). (13) 


An auxiliary scale, calibrated in 1—sin (+/2—a) is located 

immediately adjacent and parallel to the sin a@ scale. 

Once sin @ is determined, the corresponding value of 

1—sin (x/2—a) is read off on the adjacent scale and the 

sin a scale re-entered at this value. Let this point be P». 
The existing time scale may also be used, since 
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sin a 
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FicurE 3.—Proportional chart for solution of equation (11). 


sin h=cos (#/2—h)=—cos (32/2—h). 


Only the lower half of the cos / scale is used for the azimuth 
calculation. To determine the point P, at which one 





(14) 
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enters this scale, one must use either +(h—6 hr.) or 
+(h—18 hr.), whichever will give a point on the lower 
half of the time scale. 

Points P, and P, are joined by a straight line which 
intersects the uncalibrated diagonal at Ps. Then /, js 
joined with P; on the solar declination scale at the bottom 
of the nomogram and the azimuth read off at Ps on the 
azimuth scale. 

In the Northern Hemisphere, the solar azimuth is 
measured from north and, during the time of year when 
5>0, confusion may exist as to whether, for example, 
70° or 110° should be read off the azimuth scale. When 
such confusion exists, the azimuth should be determined 
for a time a little closer to noon. If the new azimuth is 
closer to the 90° point, then the correct azimuth was 70°, 
If the new azimuth is slightly farther from the 90° point, 
then 110° is the correct azimuth. 
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FORECASTING MINIMUM TEMPERATURES ON CLEAR WINTER NIGHTS 


IN AN ARID REGION 


A Comparison of Several Climatological Aids 


PAUL C, KANGIESER 


U. S. Weather Bureau Airport Station, Phoenix, Ariz. 
[Manuscript received October 7, 1957; revised July 16, 1958] 


ABSTRACT 


Minimum temperature formulas for clear nights in December, January, and February are developed by the Young 
method and tested on both original and test data. The results of these tests lead to a discussion of some weaknesses 
of the “method of arbitrary corrections.”” Jacobs’ graphical adaptation of Brunt’s equation is tested (a) using ex- 
perimentally determined values for the local soil constants to compute the “effective” soil factor, and (b) using a soil 
factor determined empirically for local meteorological data. With Brunt’s equation as a model, the physical justifica- 
tion for Young’s method is discussed, and a more direct approach suggested using the evening dry bulb and wet bulb 
temperatures and a modern method of data analysis. There is further discussion of some implicit assumptions in the 
Young method of analysis and an attempt is made to see if these assumptions are satisfied by the analysis performed 
in combining the evening dry bulb and wet bulb temperatures and the expected morning minimum using modern 
methods. In the Appendix, the application of the latter method is extended to cloudy nights and performance 
comparisons with official forecasts are made and discussed, 


1. INTRODUCTION 


Meteorologists in many parts of the world have con- 
cerned themselves during the past half century with the 
problem of forecasting the minimum temperature using 
hygrometric data [13]. Most of the methods apply best 
on clear nights, and attempt to provide the forecaster 
with an estimate of the fall of temperature to be expected 
during the night (under average radiation conditions), 
using surface temperature and hygrometric readings near 
the time of sunset. Most of the methods do not provide 
a completely objective forecast, but encourage the fore- 
caster to modify the formula value subjectively on the 
basis of expected wind and sky conditions during the 
night; hence, they may be called “climatological aids’’ to 
forecasting. 

In the Salt River Valley of Arizona, minimum temper- 
atures are of major concern to citrus growers during the 
winter season. Many of the conditions required for the 
successful functioning of most of the formulas are satisfied 
locally on the majority of winter nights because (1) the 


TABLE 1.—Data periods used in the study. 























Original Data Test Data 
Month Years included No. || Month Years included No. 
| cases || cases 
—_ = ——~ | — 
Dee __.| 1948, °50, °52, 54, ’56.| 98 || Dee.....--- 1949, "51, "53, "55___- 88 
Jan .| 1949, 51, 53, °55, °57_| 85 || Jam..........| 1950, °52, 54, ’56___-. ae 
Feb ---| 1949, ’51, °53, "55, °57.| 88 || Feb.........| 1950, °52, 54, °56.....) 71 
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skies are usually clear, (2) surface winds are light, and 
(3) the air is relatively free from pollutants. Hence, this 
should be a good place to compare the performance of 
hygrometric relationships. 

Many minimum temperature studies in the past have 
been concerned only with temperatures in the lower 
ranges, since those are the cases of primary interest in 
frost forecasting. The present study is made with the 
forecasting program of a First Order Weather Bureau 
Station in mind, and purposely includes temperatures in 
all ranges. Daily forecasts for such a station are given 
press, radio, and television distribution, so that minimum 
temperature forecasts in higher ranges, too, are of con- 
siderable operational importance. 


2. DATA USED FOR THE STUDY 


Observations made at the Weather Bureau Airport 
Station in Phoenix during December, January, and Feb- 
ruary over a 9-year period, were used in making the com- 
parisons. These data were divided into two groups as 
indicated in table 1. 

Alternate years were used in each group to reduce the 
effects of any year-to-year persistence that might be pres- 
ent in the data. Only data for clear nights were used: a 
night was called ‘clear’ if the mean sky coverage for the 
nine hourly observations (0030-0830 msT) was in the 
range 0-3/10. Hence, if any of the derived relationships 
are applied to operational forecasting, the forecaster is 
left to make the decision of whether or not the night will 
be clear. 
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Fiaure 1.—Graph of (y—d) against A for the December original data. 


3. COMPARISON OF THE METHODS 
YOUNG'S METHOD 


Some years ago, Ellison [3] made a rather extensive 
‘ comparison of a number of hygrometric formulas by ap- 
plying them to the same set of data. On the basis of that 
study, he concluded that the Young formula gave the 
best results, especially if modified by the “method of arbi- 
trary corrections.”’ Therefore, the first method tested in 
the present study was Young’s formula based on local 
data. 

Individual Young formulas for clear nights were con- 
structed using the original data for December, January, 
4 and February, respectively. Complete instructions for 

the development of the relationship are given in [3]; how- 
ever, a brief description of the procedure used will be given 
here. First a graph was plotted using (y—d) as the ordi- 
nate and A as the abscissa, where h is the 1730 mst rela- 
tive humidity percent, d is the 1730 mst dew point (°F.), 
, and y is the observed minimum temperature (°F.) the 
following morning. A straight line of best fit (minimizing 
the absolute error) was then fitted by inspection. Figure 
1 shows the graph for the December original data. Simi- 
lar graphs prepared for January and February are not 
reproduced here. 








h (1730 MST RELATIVE HUMIDITY, %) 
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Next the equation of the straight line was written in 
the form recommended in [3], known as the “Donnel 
formula,”’ 

h—n, 

No 


y=d— (1) 


where n,; is the h-intercept and n, is the negative of the 
reciprocal of the slope of the line. 

The final form recommended [3] for the equation is 
h—n, 


y=d——""4. Vet V, 


Ne 


(1A) 


where V, and V, are correction factors derived from the 
original dew point and relative humidity data, respeec- 
tively, as follows: 

1. Values of d and h for each day were written on cards, 
along with the date. 

2. Using the straight line formula (1), the formula 
minimum temperature was calculated, and the departure 
between this and the observed minimum temperature was 
written on each card. 

3. The cards were then arranged in the order of increas 
ing relative humidity. At intervals of approximately 10 
cards, the average departure of the formula minimum 
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TaBLE 2.—Formula for the minimum temperature, y, for December, 
and set of arbitrary corrections 


a 








yd Ot Vat Vs 











| “ae d Va 
15-18 | +3.5 | 13-19 +15 
19-20 415 | 29-24 0.5 
21-23 —05 | 25-28 +0.5 
24-26 +05 | 2230 | 0.5 
27-28 —-LO 31 +1.5 
29-31 —<a5 | ses | —0.5 
32-33 -—20 | 3 0 
34-38 -1.0 35-36 | +0.5 
39-47 -0.5 37-38 | —25 
48-67+ +25 | 39-464 | -0.5 





temperature estimate from the observed minimum tem- 
perature was calculated. This was then used as the arbi- 
trary correction to be applied to the formula over the 
range of relative humidity indicated by the group of cards. 
This correction was applied to the original formula estimate 
and a new departure was calculated. 

4. After the relative humidity corrections were deter- 
mined and applied, the cards were arranged in the order 
of increasing dew point. The process just outlined was 
repeated, and the arbitrary correction to be applied to 
the formula over a given range of dew point was com- 
puted. 

The formula for December and the set of arbitrary cor- 
rections are given in table 2. Similar formulas were de- 
rived for January and February. In applying the formu- 
las to the data it is not necessary to solve the equations 
each time, as it is simpler to prepare a table of values of 
(y—d) against A from the straight line, then add the cor- 
rective factors algebra‘cally to the tabular value. This 
procedure is especially helpful if a curvilinear relationship 
is drawn to the data as on figure 1. 

From these relationships, the minimum temperature 
was calculated on clear nights, using the original data 
and the independent data. The average absolute forecast 
error for each group is shown in table 4. One set of values 
was computed using the straight line formulas alone 
(without the correction factors Vz, and V,) and another 
set was computed using the complete relationship. The 
former are labeled “uncorrected” and the latter ‘cor- 
rected.” Notice that the application of the method of 
arbitrary corrections to the original data produced a 
marked increase in skill but that no increase resulted 
when the method was applied to test data. 

Next a curved line was fitted by inspection to the same 
data for each month (the line for the December data is 
shown in figure 1). The results obtained by testing each 


of these lines on both original and test data are given in 
table 4 under the heading, ‘““Young’s method, curved lines, 
Again, the method of arbitrary corrections 
Was applied to forecasts made from these lines, but this 
time the corrections were rounded off to the nearest whole 
The correction values for the December data 
Again no significant increase 


corrected.” 


deg ee, 


only are shown in table 3. 
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TABLE 3.—Correction values for Young's method, curved lines, 
December data 

















h Va d | Va 
15-19 0 13-19 | +1 
20-22 —2 20-24 0 
23-24 +1 25-28 | 0 
25-26 +1 29-30 0 
27-28 +1 31 +1 
29-31 +1 32-33 | 0 
32-33 0 34 0 
34-38 +1 35-36 || Ol +1 
39-47 0 *) Se 
48-67+ 0 39-46 | =f 





occurred when the corrections were applied to test data 
(see table 4 under heading, ‘““Young’s method, curved 
lines, corrected’’). 

To understand this lack of improvement from the use 
of the method of arbitrary corrections, one must under- 
stand what the method does. As pointed out in [8], 
application of the method to all possible combinations of 
values of d and A, using finite intervals of these, yields a 
series of straight lines that are always parallel to the 
original straight line. In [8] adjacent end-points of the 
intervals were connected with straight lines, thereby 
giving a continuous, irregular line. Such a procedure has 
no meaning, however, because there is no way equation 
(1A) can generate solutions for line segments with slopes 
different from that of the straight line given by the part 
of the formula without the correction factors. Thus a 
series of disconnected straight lines should be drawn. 

Instead of lines for all possible combinations of V, and 
V,, only the top and bottom of the series of lines were 
entered in figure 1, so that the areas within the rectangles 
contain all solutions possible with the December formula 
and tabulated corrections. The top of the farthest 
rectangle to the left (corresponding to 15 >h >18) is 
obtained from the December formula using V,=3.5, 
V.=1.5; while the bottom is obtained using the lowest 
tabulated value for Vz (—2.5). The rectangles for the 
other intervals of A and d used in the table were drawn on 
figure 1 in the same fashion. 

Notice that the general pattern of the rectangular areas 
indicates a curvilinear relationship that follows somewhat 
the curved line fitted by inspection, but there is consider- 
able up and down fluctuation which follows rather faith- 
fully the sampling variation in the original data. This 
same sampling variation is then imposed on the test 
sample, which may explain the fact that no increase in 
skill was obtained over the linear relationship; applying 
this original sampling variation to the test sample nearly 
counterbalances the increase in skill contributed by the 
general curvilinear pattern followed by the areas. Ob- 
viously, correction for the sampling variation can cause a 
significant increase in skill only if the method is tested on 
the original data alone as was done in [3]. 

Smith [10] suggested that a curve of the second degree, 
drawn on the graph of (y—d) against h, will probably 
yield better results, in the long run, than the Young 
formula. This seems borne out by the verification scores 
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TABLE 4.—Comparison of average absolute error (° F.) obtained on samples of original and test data using the various methods 





—— 





Young’s method (straight lines) 


Young’s method (curved lines) ! | 
|| Jacob’s method |} Wet bulb vs. dry 

















|} bulb method 








Uncorrected Corrected Uncorrected Corrected | | 
Original | Test | Original | Test Original Test Original Test Original Test | Original | Test 
data | data data | data data | data data | data data data data | data 
= =a - St) ae , ra er ee ee wi = See pAOBEN § allel -appnanliiaee 

| | | _ 
December 2.7 3.3 2.4 | 3.2 2.7 2.9 2.5 | 2.8 3.2 3.7 | 2.3 27 
January 3. 6 | 3.5 2.7 | 4.0 3.1 2.9 2.5 | 3.0 4.4 3.0 2.2 27 
February 3.3 3.5 | 2.5 | 3.0 2.5 | 2.8 | 2.1 2.7 3. 8 3.9 2.3 | 22 
Average 3.2 3.4 | 2.5 3.4 2.8 | 2.8 2.3 2.8 3.8 3.6 |) 2.3 | 25 





in table 1, because the curved line without correction 
factors gives a higher score on test data than the linear 
relationship with the correction factors. Note also that 
there is little difference between the scores on original 
and test data for the curved line, and on the test data for 
the curve with correction factors applied; however, the 
correction factors applied to the original data improved the 
skill considerably. This indicates that the correction 
factors are primarily correcting for the sampling variation 
in the original data, and that the curve is doing as good 
a job by itself of estimating the relationship as the curve 
with the factors applied. 


JACOBS’ METHOD 
Another aid to forecasting the minimum temperature, 
used with some success in the past, is Jacob’s [6] diagram, 
which is based on a formula for the temperature change 
at a key station in El Centro, Calif.: 


At, = 12.32(0.56—0.08 ye) yt (2) 


where ¢ is the vapor pressure (mb.), At, is the temperature 
drop (° C.), and ¢ is the time (hours) after sunset (or 
length of night). This formula was derived, in turn, 
from a formula by Brunt [1] 


9 "4 be _ 1 
201 i 1—a—bye vt. (3) 


y=T,— = 
PiCry ky 


Vr 


In (3), p; is the soil density, ¢, is the specific heat of the 
soil, and k, is the thermal conductivity of the soil, o is 
Stefan’s constant, 7’ is the radiative temperature (mean 
temperature during the night, for practical purposes), 
T, is the temperature at sunset, and y is the forecast 
minimum temperature.’ Jacobs [6] constructed a set of 
graphs based on formula (2) and provided factors to 
correct for variations in 7’ and for different values of the 


? 





“effective” soil factor (s), where s=yxp,er/k;. a and 6 are 
constants in an empirical equation developed by Brunt [1] 
which expresses the relationship between the surface 
vapor pressure, the total long-wave radiation downward 
from the atmosphere, and the total black-body radiation 
at a given surface temperature on clear nights. The 
values of a and 6 used here are those used by Jacobs 
(a 0.44, 6 = 0.08). 


1 Jacobs used 7 for the minimum temperature, but y will be used here to be consistent 
with the notation used above. 





Local average values of p;=2.70 gm. em.~* and ¢,;=0.20 
cal. gm.~' °C.~' were obtained from the College of Agri- 
culture at Arizona State University. The mean value of 
k,=4.7X107* cal. sec.~' em.~' °C.-' given by Jacobs [5] 
was used, along with a value of t=14.0 hr. The shelter 
temperature at 1730 mst was substituted for 7, and the 
vapor pressure at the same time was used for e. 

In computing formula (2), Jacobs used a value of 
T=280° AA. In the present study, a value of 7=283° 
AA was used since the mean nighttime temperature at 
Phoenix during the winter averages near 50° F. Actually, 
the numerical value of 7 used in developing the basic 
formula is arbitrary, since a correction is ultimately in- 
cluded for the observed (or forecast) mean temperature 
during the night; it becomes important only if the corree- 
tion is not applied. 

The values of the variables and constants given above, 
gave the following formula for Phoenix: 


T,—y=16.1(0.56—0.08y@) yi. (4) 
This equation expresses (7',—y) in degrees Celsius. Curve 
A on figure 2 gives (7,—y) in degrees Fahrenheit, using a 
constant value of t=14.0 hours. Figure 3 gives the value 
of the correction factor (f) by which the value of (7,—y) 
from figure 2 must be multiplied to correct for the actual 
mean temperature during the night. In the verification, 
the observed value of T=%(7,+y) was used as the mean 
temperature during the night. Formula (4) was verified 
only on the December original data, since it became ob- 
vious during the verification that the formula was con- 
sistently forecasting minimum temperatures that were 
much too low. An average error of 11.5° F. was obtained. 
Brunt [1] states that the chief difficulty in the way of 
using formula (3) for forecasting the night minimum tem- 
perature lies in the uncertainty of the value of the coeffi- 
cient p,¢;Vyk;, and that the addition of 20 percent of water 
to dry soil increases this factor five-fold. Jacobs [6] sug- 
gests that these factors can probably best be derived em- 
pirically, so equation (3) was written in the form 





Ses (y— T,) 

2oT*yt Z 
DS —_ — ——’ (5) 
pitivky (1—a) —bye 


From the values stated above of o, 7, t, a and 6, and ob- 
served values of y, 7, and e (1730 mst) for all clear nights 
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Fievre 2.—(7,.—y) plotted against vapor pressure (e) at 1730 mst. 
Curve A=curve for equation (4). Curve B=curve for equation 
(6). 


in December, January, and February of the original data, 
a value of 1/p,c: vk, was computed for each of the 271 
nights. A frequency polygon was constructed from these 
values and the median value computed. The distribu- 
tion of values turned out to be very regular (little or no 
skewness) so that the mean and the median coincided at 
33.5. From the frequency polygon, a percentile estimate 
of the standard deviation using 10 percentiles [2] gave a 
value of 7.5. This estimate (statistical efficiency 0.92) of 
the standard deviation is included only to give some idea 
of the spread of the data about the mean value. 

The median value of 1/p,c, Vk, was then substituted into 
equation (3) and by rearrangement a new equation similar 
to (4) was obtained: 


T,—y=11.1 (0.56—0.08yé) ¥2. (6) 


Like equation (4), this equation gives (7',—y) in degrees 
Celsius, assuming a mean temperature during the night of 
283° AA. Curve B on figure 2 gives expected values of 
(T,—y) in degrees Fahrenheit, using a constant value of 
t=14.0 hr. Figure 3 may be used in the same manner as 
with equation (4) to correct for the observed mean tem- 
perature during the night. 

As in the case of equation (4), the verification of this 
equation was carried out using observed values of 7'= 
4(T,+-y). As can be seen in table 4, the results are con- 
siderably improved over those of equation (4). The use 
of Jacobs’ equation on an operational basis requires, in 
addition to a forecast of sky condition, an estimate of the 
mean temperature during the night; however, using the 
observed 1730 mst temperature and a rough estimate of 
the minimum temperature to compute 7=\(7,+-y) would 
probably have given results on the same data nearly com- 
parable to those obtained here. 


A MORE DIRECT APPROACH 


A climatological aid for forecasting the minimum 
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Figure 3.— Correction factor (f) by which (7',—y) from equations 
(4) and (6) should be multiplied to correct for the mean tempera- 
ture during the night. 


temperature should be a simple relationship based on 
evening data, which the forecaster can use for guidance in 
making his forecast. Most of the methods developed in 
past years have used evening temperature and moisture 
parameters to forecast the minimum temperature the 
following morning, and physical justification for this pro- 
cedure can be found in Brunt’s equation. In equation (3) 
all variables are assumed to be constants with the excep- 
tion of y, T,, T, t, and e. We have seen, however, that 7 
may be closely approximated by (y+ 7,)/2, and t may be 
taken as a constant (14.0 hr. was used above). Hence, 
the only variables left are y, 7,, and e. Young’s method 
essentially substitutes two other functions of 7, and e; 
namely, h and d. Another possibility for a 3-variable 
combination that seems quite feasible, in terms of Brunt’s 
equation, is y, 7,, and the wet bulb temperature 7’,. 
Both 7, and 7), are measured directly in the instrument 
shelter and so are somewhat easier to use on an operational 
basis than d and h. 

Given three variables, 7,, 7,, and y, between which a 
joint relationship is suspected, there are a number of ways 
to approximate this relationship [4], [7]. One of the most 
straightforward, however, is to plot 7, along one coordi- 
nate axis, 7, along the other, and the numerical value of 
y at the intersection point in the field of the graph corre- 
sponding to particular values of 7,, 7,. Isopleths of 
constant values of y may then be drawn by one of several 
methods [4], [7]. A diagram of this type is shown in 
figure 4, based on the 271 cases of original data for Decem- 
ber, January, and February combined. The isopleths 
were drawn by the method of substratification and suc- 
cessive graphical approximation described in [4}. 

The advantage of this type of combination of variables 
is that no initial assumption is made about the form of the 
relationship, and the data are given freedom to indicate 
their relationship within the limitations of the sample 
size. This is in contrast to the Young formula, for ex- 
ample, which plots two variables along one axis as (y—d), 
and thereby implicitly assumes that these two variables 
have equal weight in the relationship with A, and that the 
regression of y on d (or vice versa) is linear for constant A. 
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Fraure 4.—Graph showing the estimated joint relationship between 7',, T., and y on clear nights. 


Lines of constant d and h were added 


using hygrometric table after analysis was performed. 


These assumptions become readily apparent if the relation- 
ship between y, d, and h is examined on a graph using 
(y—d) as the ordinate and A as the absicssa. Such a 
relationship can be expressed in the following general 
functional form, 

y—d=f(h) (7) 
or, 


y=f(h)+d (8) 


where {(h) is assumed to be a continuous function for all 
values. 
Now consider the effect of holding h constant, 


y=constant+d. (9) 


The linearity between y and d for constant h and the 
assumption of equal weight of y and d with respect to h 
are shown by (9), and can be graphically illustrated on a 
diagram of the kind shown in figure 4, by using y as the 
ordinate, d as abscissa, and h in the field of the diagram. 
Such a diagram will always have the following set of char- 
acteristics no matter how complicated the curve drawn 


on figure 1 to indicate the relationship between y, d, and h, 
lines of constant A will always be straight lines, and the 
slope of these lines will always be 1. As an example, the 
relationship between y, d, and h shown by the curved line 
in figure 1, is converted into this kind of diagram in figures 
(solid lines). 

Since the Young method has the implicit assumptions 
mentioned above, it might be of interest to see how wel 
these assumptions are satisfied by the analysis in figure 4, 
since that analysis was carried out unencumbered by thos 
assumptions. To do this, however, we must convert the 
relationship shown in figure 4 into one between y, d, andh. 
This can be done simply by drawing lines of constant 
and A in figure 4 using a hygrometric table. These lines cal 
then be used to graph the relationship between y, d, and| 
in figure 5 (dashed lines), for a comparison with the lines 
(solid) drawn using the curvilinear relationship in figure |. 
Obviously, such a comparison can be used only W 
indicate differences between the results of the two methods 
on this sample of data, since neither one is known to be 
“correct.”” The dashed lines, which were free to assume 
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-Graph showing relationship between d, y, and h for clear 
Solid lines show the relationship represented by the curve 
Dashed lines show the relationship from the analysis 


Figure 5. 
cases. 
in figure 1. 
on figure 4. 


any configuration in figure 5, seem quite straight for values 
of y less than about 40°. Thus the assumption of linearity 
of the regression of y on d for constant h seems fulfilled 
by the figure 4 analysis for lower values of y, but seems to 
break down at higher values of y (and higher values of 
dand h). 

Young [12] found that forecasts made from a straight 
line relationship on a graph like figure 1 were considerably 
in error at high values of d and A (and also for very low 
values) and developed the method of arbitrary corrections 
to correct for this. It should be emphasized here, how- 
ever, that the statement above about the indicated lack 
of linearity in figure 5, should not be confused with 
Young’s finding. The lack of linearity indicated in 
figure 5 is between the regression of y on d for constant h, 
and if real, cannot be corrected by any method as long as 
the basic graph is like that used in figure 1. Obviously, 
the method of arbitrary corrections also operates within 
the framework of the two assumptions described by 
equation (9). 

For the dashed lines in figure 5, the slope is approx- 
imately 1.3 for values of y less than about 40. For higher 
values of y, the slope becomes variable. 


Thus the assumptions implicit in the Young method 
are only approximately fulfilled by the original analysis 
in figure 4. Although we cannot draw definite conclusions 
from the comparison in figure 5, it is possible that some of 
the difference in skill on the test data between the analysis 
in figure 1 and that in figure 4 (see table 4) is due to a 
of conformity of the data to the assumptions made 
In fgure 1, 


To provide some standard of comparison for the per- 
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TaBLE 5.—Comparison of average absolute error for the climatological 
aid (T., vs. T,) and official forecasts for clear nights 











Season No. cases | T, vs. Ty | Official 

forecasts 
1954 ‘ ee ; "reais 63 2.7 3.8 
1956 : 57 2.2 3.2 
1958 53 2.2 22 





formance of figure 4, temperature estimates made using 
1730 mst data were compared with official forecasts made 
by the staff at WBAS Phoenix. The official forecasts 
of the minimum temperature the following morning are 
made on an operational basis at 1930 msr each evening. 
The average absolute error by seasons is shown in table 
5 for all clear cases. The seasons used were from the 
test sample for figure 4 (1954 season=Dec. 1953, Jan. 
and Feb. 1954; 1956 season=Dec. 1955, Jan. and Feb. 
1956; etc.). It should be pointed out in this comparison, 
that the climatological aid “knew” that the night would 
be clear, while the staff forecasters did not. 

It is of interest to note that the climatological aid 
averaged approximately 1° F. better than the official 
forecasts in the first two seasons, but in the last season 
the skill was about the same for both. During the 1958 
season the Meteorologist in Charge placed the clima- 
tological aid in the station forecast manual and the fore- 
casters used it on a discretionary basis. There is no way 
of determining what part (if any) of the improvement in 
skill shown in the 1958 season was due to its use. 


4. CONCLUSIONS 
YOUNG'S METHOD 


In this method a graphical combination of y, d, and 
h is performed, but the combination is done within the 
framework of the following assumptions: (1) the regres- 
sion of y on d, for constant A, is linear, and (2) y and d 
have equal weight in the relationship with A. If these 
assumptions are not fulfilled by the “population,” a loss 
in skill can be expected over that obtained by a method 
of analysis allowing greater degrees of freedom. 

Another logical combination of variables, in terms 
of Brunt’s equation, is 7,, 7T,, and y. An advantage 
of this combination is that the independent variables 
(T, and T,) are measured directly in the instrument 
shelter. These variables were combined on a graph 
using 7, as ordinate and 7, as abscissa. If estimates 
of y are desired in terms of d and h, as in Young’s method, 
it is only necessary to add lines of constant d and h to 
this diagram with a hygrometric table. 

With Brunt’s equation as a working physical model, 
Young’s method which combines A, d, and y is seen to 
be essentially the same as the method employing 7, 7'y, 
and y. It is possible that the improvement in skill shown 
on test data by the latter method is largely due to the 
use of an improved method of data analysis. Whether 
or not this is true, we have seen that the Young method 








) 


of combining d, A, and y introduces assumptions for 
which no a priori physical reasons have been given. 
There seems little reason for imposing these assumptions 
today, when we have a number of good methods for 
estimating the joint relationship between several variables 
which do not require such arbitrary restrictions. 

A relationship between y, d, and h was derived using a 
method of analysis that did not arbitrarily place these two 
restrictions on the data. The following conclusions were 
drawn for the particular analysis performed on this par- 
ticular set of data: (1) assumption No. 1 was well satisfied 
for observed minima below about 40° F. but broke down 
at higher minima, and (2) assumption No. 2 was not quite 
as well satisfied, even at temperatures in the lower ranges. 

If the Young method of combining y, d, and h is used, 
it is recommended that instead of preserving the original 
first-degree estimate of the plotted relationship between 
(y—d) and h by applying the corrections in the form of 
formula (1A), a curvilinear relationship of the second 
degree be drawn. Forecasts made from this line will 
probably be subject to less loss in skill caused by imposing 
original-data sampling variation on future samples. 


JACOBS’ METHOD 


This method does not seem to be more capable of uni- 
versal application (application without determining local 
constants by data processing) than other methods. This 
may be due primarily to the difficulty of obtaining reliable 
estimates of the soil factors, which probably not only 
undergo considerable seasonal variation, but also may 
differ, on the average, between clear and cloudy nights. 


A MORE DIRECT APPROACH 


Inspection of Brunt’s equation (3) provides a physical 
justification for using the evening dry bulb temperature 
(T,) and the evening vapor pressure (e) to forecast the 
minimum temperature the next morning (y). Since dew 
point (d) and relative humidity (h) are functions of T, 
and e, the combination of d, h, and y also seems logical, 
but it is recommended that a method of estimating the 
joint relationship between these variables be used that 
does not place unnecessary restrictions on the data. 


5S. SUGGESTIONS FOR FURTHER RESEARCH 


This study has discussed methods which, in effect, use 
T, and e to estimate y. Obviously, there are many other 
variables of great importance, but the introduction of 
these into the relationship was beyond the scope of this 
project. Some of these important additional variables are 
(1) wind, (2) turbulence, (3) advection, (4) soil moisture, 
(5) soil temperature, (6) air temperature aloft, and (7) 
moisture distribution aloft. One approach to future re- 
search in this field might be to incorporate these variables 
into the relationship, perhaps using the “residual method”’ 
9], which has the advantage of showing whether an addi- 
tional variable improves a forecast based on other vari- 
ables. 
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Another approach may be to work directly with Brunt 
expression for the fall in temperature after sunset, which js 


2 fy 


vr pier ky ’ 

where Fy is net radiation, and ¢ is the time after sunget 
Since changes in soil moisture are mainly responsible for 
changes in the factor p,¢;yk, in a given locality, it may be 
possible to substitute a soil moisture reading, empirically, 
for this factor. Ry can be measured using fairly inex. 
pensive equipment [11]. One procedure would be to take 
a reading of Ry near or shortly after sunset and, following 
Brunt’s suggestion, consider it to be constant throughout 
a given clear night. Forecasts made in this fashion woul 
have to be modified subjectively for advection, wind, and 
turbulence, just as those discussed previously. 
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APPENDIX 
EXTENSION OF THE METHOD IN FIGURE 4 TO CLOUDY NIGHTS 


To extend the usefulness of the climatological aid, 
graphs similar to figure 4 were prepared for two other 
classes of nights, as follows: 

Class 1. This class contains principally nights, with 
high broken to overcast conditions and with no persistent 
ceiling below 10,000 feet. Criteria: on the morning § 
following the evening on which the forecast was made, 
the observed mean sky coverage on the nine hourly 
observations 0030-0830 mst, inclusive, was in the range 
4/10—10/10 and a ceiling of 10,000 feet or less was reported 
on no more than 4 of these observations. 

Class 2. This class contains principally nights with 
broken to overcast skies and with a persistent ceiling 
below 10,000 feet. Criteria: on the morning following 
the evening on which the forecast was made, the observed 
mean sky coverage on the nine hourly observations 
0030-0830 mst, inclusive, was in the range 4/10—10/I0 
and a ceiling of 10,000 feet or less was reported on 4 0 
more of these observations. 

It will be noted that on figure 4 the minimum temper 
ature lines are not extended into the region of the graph 
above a relative humidity of 90 percent. Extension of 
the lines into that region would have little meaning 
because a relative humidity that high in the early evening 
usually leads to ground fog by morning, automatically 
placing the case in the “cloudy, class 2” category. For 
this same reason there were practically no observed case 
in this region of the diagram for the “cloudy, class 1’ 
category, (fig. A-1) so the minimum temperature line 
above a relative humidity of 90 percent were omitted 
this diagram also. 
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Ficure A-1.—Graph showing the estimated joint relationship 
between 7, 7, and y on cloudy nights (class 1). Lines of con- 
stant y (estimated) are drawn in the field of the diagram. 


On the “cloudy, class 2” diagram (fig. A-2), the 
minimum temperature lines were drawn from the inter- 
section with the 90 percent relative humidity line to a 
point of intersection with the 100 percent relative humidity 
line, corresponding to a dry bulb temperature about 2° 
higher than the expected minimum temperature line. 
This was an arbitrary decision, based on the assumption 
that if the relative humidity were 100 percent at 1730 
ust, the dry bulb temperature would only fall about 2° 
during the night. This decision was based on a few 
observed cases on foggy nights, but the frequency of such 
nights in Phoenix is very low so the evidence for deter- 
mining the configuration of the lines in this region of the 
chart was flimsy. Physical reasoning would indicate, 
however, that as the relative humidity approaches 100 
percent the release of the latent heat of condensation, as 
fog formation begins and progresses, should reduce the 
rate of temperature drop materially, and the minimum 
temperature lines at high relative humidities should be 
drawn with that fact in mind. Enough data might be 
available at coastal stations to do a much better job of 
analysis in that region of the chart. 

Because of the fewer number of cases in the two cloudy 
categories, the original data sample was extended to 
include the following years: December 1948-1955; January 
and February 1949-1956. This gave 175 cases in class 1, 
and 85 cases in class 2. The verification on original 
data is shown in table A-1. 

The verification on test data was carried out for 


December 1956, 1957; January and February 1957, 1958 
(a total of 64 cases for class 1 and 21 cases for class 2). 
Because of the scarcity of cases, the verification in table 
A-1 is given only for the whole test sample for both the 
climatological aid and the official forecasts. 

The importance of being able to forecast cloudiness 
the night ‘s especially apparent in the great 
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Figure A-%.—Graph showing estimated joint relationship between 
T., Tw, and y for cloudy nights (class 2). Lines of constant y 
(estimated) are drawn in the field of the diagram. 


differences in skill between the climatological aid and the 
official forecasts for class 2, since the former essentially 
“knew” the ensuing sky conditions at the time of 
“forecast.” 


TaBLE A-—1.—Average absolute error for cloudy cases using original 
data and comparison with official forecasts using test data 





Original data Test data 





| Class 1 Class 2 Class 1 Class 2 
Sak < Sees eee Si ot ae _ar <5 ten : 
| used | 
| No. | Avg.| No. | Avg.|| No. T. | Off. | No. -. oft. 
| cases | abs. | cases | abs. | cases vs. | fore- | cases! vs. | fore- 
error error Tw | casts Tw | casts 
Dec 48-55 5 3.0 28 23 64 3. 0 3.7 21 3. 6 6.4 
ae... 49-56 65 | 3.0 36] 25 
Feb .| 49-56 a} oes 2.6 
Average 175 3.0 8&5 25 
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SOLAR ENERGY AND SUNSHINE HOURS AT ATHENS, GREECE 
GEORGE JAC. MACRIS 


National Observatory of Athens, Greece 


{Manuscript received November 3, 1958; revised December 15, 1958] 


1. INTRODUCTION 


Measurements of solar radiation at the earth’s surface 
are restricted to widely scattered locations and in many 
cases represent relatively short periods of record. Meas- 
urements of the hours of sunshine, on the other hand, 
are more extensive in both time and space. 

A number of investigators have successfully applied 
various forms of the linear relationship between the total 
daily amount of solar radiation incident on a horizontal 
surface and the hours of sunshine as proposed by Angstrém 
(|. This relationship may be expressed as follows: 

: 
@ =a+b S 
Yo Do 
where Q=total amount of solar radiation impinging on 
a horizontal square centimeter at the earth’s 
surface per day (expressed in ly. per day). 
Q=maximum total amount of solar radiation 
impinging on a horizontal square centimeter at 
the earth’s surface on a perfectly clear day (ly. 
per day). 
S=number of hours of sunshine measured on a 
sunshine recorder. 
Sy=maximum number of hours of sunshine re- 
corded under perfectly clear conditions. 
a and 6 are constants. 

The application of such a formula makes it possible 
tomake a quantitative estimate of solar radiation amounts 
in cases where only the more common sunshine measure- 
ments have been made. For example, at Athens, Greece 
the record of sunshine hours extends from 1901 to the 
present, while solar radiation measurements have been 
made only during the last few years! [i2]. Thus the 
linear relationship cited above can be used to extend the 
estimate of solar radiation and its frequency distribution 
over a far longer period than actual records cover. 


2. INSTRUMENTS 


At the National Observatory, Athens, Greece [12] 
daily totals of solar radiation are measured by a Gorc- 
azinsky pyrheliometer (No. 111, Richard). For compari- 


_ 


Some carlier investigations of solar energy at Athens, Greece were reported by Kara- 
plperis 


son a Kipp and Zonen Solarimeter (No. 604) is used. 
The number of hours of sunshine is recorded by a Camp- 
bell-Stokes Sunshine Recorder. Four years of simulta- 
neous measurements of Q and S, 1953 through 1956, are 
included in this study. 


3. CALCULATION OF THE LINEAR RELATIONSHIP 


Y 


; S 
Daily values of @ versus — were plotted for each month 
Ad) 0 


and the linear equation for each month was determined 
by the method of least squares. Examples of such 
scatter diagrams for the months of March and August 
are shown in figures 1 and 2. 

A lack of scatter in the data for the months of June, 
July, and August, due to persistent clear conditions 
made the resulting linear relationships for these months 
of doubtful value. This difficulty was overcome by 
approximating the monthly values of 6 by the average 
value of 6 for the other nine months. 


4. RESULTS 


The monthly values of the constants a and 6 determined 
in the above manner are given in table 1. These values 
give an average of a=.34 and b=.63. These compare 
with the values found by Fritz and MacDonald [4] for the 
United States of a=.35 and 6=.61. 

The variation from month to month of a and 6 and the 
variation from location to location in annual values of 


TABLE 1.—Monthly values of solar energy (Qo and Q/Qy), sunshine 
(So and S/So), and linear regression coefiicients (a and b) of Q/Q, vs 
S/So, Athens, Greece, 1953-56. Qo is in ly./day, So in hours/day 


Energy Sunshine Inter- Linear 
Month : 4 cept Slope 
a b 
Qo Q/Q Si S/S 

January 283 65 9.6 43 . 39 t) 
February 365 72 10. 4 49 39 67 
March 559 61 11.5 fT) 2 65 
April 565 83 12.3 67 37 68 
May 628 83 13. 2 77 34 63 
June 701 s4 14.2 87 2 63 
July 689 SY 14.2 W6 2 63 
August 6565 &5 13.0 6 25 63 
September 5M 78 12.5 80) 32 as 
October 436 67 10.9 62 30 i) 
November 277 74 9.4 52 41 (4 
December 2265 76 sf 55 44 58 
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Figure 1.—Scatter diagram and linear regression line of Q/Qo vs S/So for March, 1953-1956, Athens, Greece. 
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Ficure 2.—Scatter diagram and linear regression line of Q/Qo vs S/S» for August, 1953-1956, Athens, Greece. 
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Figure 3.—Annual march of solar insolation at top of atmosphere 
at latitude of Athens, Greece (curve 1). Annual march of maxi- 
mum solar insolation possible at surface on clear days at latitude 
of Athens, Greece (curve 2). 


these constants are discussed by Angstrom [2], Fritz and 
MacDonald [4], Fritz [5], Kimball [8], [9] and others [6], 
(7). Of particular importance in this case are marked 
seasonal variations in the physical characteristics and 
distribution of clouds. 

Table 1 also contains monthly values of S, and Q, as 
determined from the examination of the maximum 
perfectly clear day values of S and Q recorded during each 
month during the period of record. 


5. DEPLETION OF SOLAR ENERGY BY THE 
ATMOSPHERE 


Curve 1 of figure 3 is the plot of average daily totals of 
solar radiation outside the atmosphere at the latitude 
of Athens, from data by Bernhardt and Philipps [3]. This 
may be compared with curve 2, a plot of S,, to illustrate 
the monthly variation in the depletion of solar radiation 
by a “clear” atmosphere. The atmosphere is clearer and 
also drier during the fall and winter than in the summer. 

During the months of April and May, the radiation 
at the surface is lower than expected as curve 2 indicates. 
This may be due to the additional absorption of solar 
energy by water vapor in the moist air masses which move 
off the Mediterranean from the southwest, and which are 
prevalent in these two months. In other months of the 
year the prevailing winds are northeasterly [11] and the 
water vapor content of the atmosphere is much lower. 


6. FREQUENCY DISTRIBUTIONS OF 
TOTAL HOURS OF SUNSHINE 


Frejueney distribution histograms of the total hours 
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Figure 4.—Histogram of total hours of sunshine for January, 
1953-1956, Athens, Greece. 
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Figure 5.—Histogram of total hours of sunshine for July, 1953- 
1956, Athens, Greece. 


of sunshine were determined for Athens on the basis of 
records from 1901 until 1956. Two examples of these 
histograms from 1953 to 1956 are illustrated in figures 4 
and 5. Of particular importance is the frequent occur- 
rence of clear days during the summer months. In July 
87 percent of the days have more than 13 hours of sun- 
shine. This is in sharp contrast with the more nearly 
uniform distribution of the winter months with January 
showing 20 percent of its days with less than one hour of 
sunshine. 
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THE WEATHER AND CIRCULATION OF JANUARY 1959 
A MONTH OF EXCEPTIONAL PERSISTENCE FROM THE PRECEDING DECEMBER 


L. P. STARK 


Extended Forecast Section, U.S. Weather Bureau, Washington, D.C. 


1. PERSISTENCE 
FROM DECEMBER 1958 TO JANUARY 1959 


WEATHER ANOMALIES 


Very high persistence over North America from Decem- 
ber to January was one of the outstanding facets of the 
monthly mean weather and circulation of January 1959. 
One should anticipate a normal increase in month-to- 
month persistence of temperature, mean 700-mb. heights, 
and to a lesser degree, precipitation in passing from fall 
into winter [1]. However, stability of pattern and weather 
anomalies from December to January was so great this 
year that special emphasis is warranted. 

Class changes of monthly mean temperature in the 
United States are shown in table 1. These data were 
derived from an array of 100 well-distributed stations and 
represent changes in classes based on the normal distribu- 
tion of temperature at each city. The classes much above 
and much below normally occur 12! percent of the time 
each, and above, normal, and below normally occur 25 
percent of the time each. On examining the table, one 
is first struck by the absence of any change exceeding two 
But even more significant, and certainly more 
indicative of the month’s persistence, is the fact that 96 
percent of the stations did not change by more than one 
class. 


classes. 


It is interesting to compare the results in table 1 with 
the averages determined by Namias [1] for the period 
1942-1954. In his study the average number of 0+1 


TaBLe 1.—Class changes of weather anomalies in the United States 


from December 1958 to January 1959 
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class changes (no change or a change of one class) for 
December to January was 71 percent. By comparison, 
January 1959 with 96 percent was phenomenally static. 
Fairly high persistence has been found in other Januarys 
(e.g., 1956 [2]), but in the last 17 years there has not been 
a January with smaller temperature changes than this 
year. Parenthetically, it should be noted that 0+ 1 class 
temperature persistence greater than 90 percent has been 
observed only 17 times since April 1942 (8 percent of the 
time). 

Precipitation from December 1958 to January 1959 was 
also quite persistent. Table 1 shows a 47 percent value 
of the 0 class change in precipitation for January. This 
represents rather high stability when compared to the 
December-January average of 37 percent for the years 
1942-1954 [1]. In fact, there were only three previous 
Januarys in the past 17 years in which the 0 class change 
of precipitation was greater than the 47 percent of 1959. 

It is difficult to make direct comparison of persistence of 
precipitation with that of temperature since values are 
computed for 0 class change for the former and 0+-1 class 
change for the latter. However, the two can be compared 
by expressing them in terms of their long-period average. 
For example, this month’s persistence of precipitation was 
127 percent of the 13-year average value, while the cor- 
responding figure for temperature was 135 percent. 


CIRCULATION 


Persistence of the 700-mb. pattern was as pronounced 
as that of the weather anomalies. A qualitative com- 
parison of the 30-day mean 700-mb. chart for January 
1959 (fig. 1) with that of December 1958 (fig. 1 of [3)]) 
shows a ridge in western North America, a trough along 
the east coast, and a Low in Canada. Height departures 
from normal in the same figures show a remarkable 
similarity. 

A quantitative evaluation of circulation persistence 
can be expressed in terms of a lag correlation between the 
patterns of height anomalies. This was done earlier [1] 
for the period 1933-1950, using a grid extending from 30° 
to 50° N. and from 70° to 130° W. An average lag cor- 


relation of +0.36 was obtained for all December to Janu- 
ary periods. 


A similar correlation computed for Decem- 
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Ficure 1.—Mean 700-mb. contours (solid) and height departures from normal (dotted) (both in tens of feet) for January 1959. Heavy 
dashed arrows represent the axis of the mean 700-mb. wind maximum. Widespread positive height departures over the polar region 


reflected intense blocking there. 


ber 1958 to January 1959 was +0.61, substantially greater 
than the average for the 18-year period above. 

There was one feature of the circulation not wholly 
consistent with the high persistence. A marked change 
occurred in the strength of the 700-mb. temperate wester- 
lies (35°-55° N., computed twice monthly for overlapping 
30-day periods). The zonal index for the western half of 
the Northern Hemisphere for December was 12 m.p.s., 
0.7 m.p.s. above normal. During January the index de- 
creased to 9.9 m.p.s., 1.9 m.p.s. below normal. This is the 
first month in the last nine in which the temperate index 
averaged below normal. The value of the 30-day mean 


sea level index in temperate latitudes averaged 1.6 m.p5. 
below normal, following a near normal index in December. 
The configuration of systems typical of low index can be 
seen on the 30-day mean sea level chart (fig. 2). These 
features include split Aleutian and Icelandic Lows, strong 
polar Highs, and weak subtropical Highs. 


2. MAJOR FEATURES 
OF THE NORTHERN HEMISPHERE CIRCULATION 


The hemispheric circulation in January 1959 (fig. }) 
consisted of a wave number of four in middle latitudes. 
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The t: ughs at higher latitudes in the western half of the 
North: rn Hemisphere were rather ‘weak compared with 
those 0! the eastern half. This is shown by the relative 
intensi'v of the 700-mb. 30-day mean height departures 
from normal. The troughs near North America were less 
intense and did not penetrate as far northward as did 
those of Europe and Asia. 

The circulation in midtroposphere in far northern lati- 
tudes of the Western Hemisphere was much more anti- 
eyelonic than that which normally prevails [4]. The 
large High over the Beaufort Sea was an especially sig- 
nificant center of action in influencing weather in the 
United States. It was part of a very extensive band of 
positive height anomalies which blanketed the area from 
the Bering Sea to Denmark Strait with heights ranging 
from 400 to 700 feet above normal. The unusual circula- 
tion in that region, which included anomalous centers of 
the greatest absolute magnitude in figure 1 and which 
dominated the flow in the Northern Hemisphere, was 
also conspicuous in the weather and circulation regime of 
December. In fact, high-latitude blocking over North 
America has been a circulation feature of varying magni- 
tude in each month since June 1958. 

A fairly well-developed Low in the Gulf of Alaska and 
its attendant trough south-southwestward maintained 
the approximate intensity observed in December (fig. 1 of 
(3}), but important changes occurred. ‘The trough ac- 
quired a positive tilt as the lower portion retrograded 
some 25° from its December position. As the trough 
along the coast of Asia deepened, the wavelength in the 
Pacific shortened. The fast westerlies observed in mid- 
Pacific in December were replaced by a difluent zone near 
the International Date Line, and the mean 700-mb. wind 
speed maximum (shown by a dashed arrow in fig. 1) was 
displaced 10° southward in January from its position near 
37° N. in December (see fig. 3A of [3]). 

Though the ridge-trough couplet over the United States 
and Canada appears to have been of small amplitude, the 
flow was vigorous enough to keep the West warm and the 
rest of the United States quite cold (fig. 3). The 30-day 
mean 700—1000-mb. thickness anomaly field (fig. 4) shows 
the extent of the cold air and its extremely cold pool over 
northern British Columbia. The channel of cold air from 
northwestern Canada to Tennessee was practically coin- 
cident with a strong ridge of high pressure at sea level 
(fig. 2) and with the path of maximum anticyclonic activ- 
ity, illustrated in figure 2 as a heavy dashed arrow. (See 
also tracks of centers of anticyclones at sea level, Chart IX, 
of [5].) 

The ridge over Greenland, principally a development of 
the first half of the month, was closely associated with the 
strength of the trough in Europe. The westerly flow into 


central Europe met the northwesterly flow from the Arctic 
ina zone of strong confluence, and cold air was thus con- 
tained north of 50° N. in the mean (see fig. 4). 

Stroig westerlies across northern Asia, with the axis of 
the mean 700-mb. wind speed maximum reaching 60° N. 
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Figure 2.— Mean sea level isobars (in millibars) for January 1959. 


Heavy dashed arrows represent the main anticyclone track in 
January. The prominent ridge from the Arctic Ocean to the Gulf 
of Mexico shows the southward drive of cold air. 
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Ficure 3. 


Departure of average temperature from normal (° F.), 
January 1959. Hatching indicates areas of normal or above 
normal temperatures. The barrier effect of the Rockies is well 
shown as cold air stayed east of Divide. 


(fig. 1), accommodated the long wavelength downstream 
to the next trough in eastern Asia. 
supported the strong, full-latitude trough observed in that 
area, which is also a zone of high trough frequency clima- 
tologically [6]. 


Adequate cold air 


3. VARIABILITY IN THE WESTERN HEMISPHERE 
DURING JANUARY 


Since the 30-day mean smooths shorter-term variations 


in the pattern, a breakdown of the mean into 15-day com- 
ponents is of some interest. 
the January pattern closely resembled that of the preced- 
ing month. 
sistence within the month. 
circulation of the first and last halves of the month (figs. 5 
and 6) reveals rather sudden changes. 


It has been noted above that 


One might thereby reasonably expect per- 
But a comparison of the 
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Ficure 4.— Mean thickness (1000-700 mb.) departures from normal 
(tens of feet) for January 1959. Below normal thickness areas 
(shaded) were almost identical to areas of below normal surface 
temperature (fig. 3). 


A downstream transfer of westerly momentum from the 
coastal trough near Asia assisted in the erosion of the sub- 
tropical ridge in the mid-Pacific, where heights fell as much 
as 1000 feet (fig. 7). This allowed retrogression of the 
trough in the Gulf of Alaska (fig. 6), where heights in- 
creased over 300 feet, and its consolidation with the Kona 
trough west of Hawaii. 

The moderately long wavelength between the trough in 
the eastern Pacific and that in the western Atlantic was 
supported during the first haif of January by rather strong 
westerlies. But the half-wavelength downstream from the 
ridge along the west coast was considerably longer than 
that usually observed, a condition which usually precedes 
a circulation change of major proportions. The vigorous 
new trough over the Mississippi Valley in the succeeding 
half-month (fig. 6) is an excellent example of a new mean 
trough development as described by Namias [7] and ap- 
parently was associated with: 

a. A wavelength too long to be supported by the speed 

of the westerlies. 

b. Retrogression of the ridge in western North America 

and its intensification as it amalgamated with the Arctic 

Basin High. 

c. Strengthening of the trough in central Canada. 

d. Northeastward extension of the trough in Mexico. 

This classical development of a new trough is also a 
prime example of amplification of the contour pattern. 
The intensity of the change from figure 5 to figure 6 is 
depicted in figure 7, the change in height anomaly from 
the first to the last half of January. Meanwhile the trough 
in the western Atlantic was swept eastward to the Azores, 
where heights decreased over 400 feet. South of New- 
foundiand height increases exceeded 700 feet as a strong 
ridge replaced the trough. 

Temperature anomalies for each half of the month, as 
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Friaure 5.—Mean 700-mb. contours (solid) and height departures 
from normal (dotted) (both in tens of feet) for January 1-1}, 
1959. Fast westerlies and a long half wavelength dominated the 
United States. 


seen in figure 8 A and B, appear very much alike in spite 
of differences in circulation. However, an examination 
of figure 8C shows significant changes related to the evoly- 
ing circulation. Note especially the negative anomalous 
temperature change in figure 8C and compare it with the 
height change in figure 7. Positive height anomalies of 
January 1-15 were replaced by negative anomalies of 
January 16-30 from Hudson Bay to the Gulf of Mexico 
and from the Continental Divide to the Ohio Valley. 
At the same time, temperature anomalies became more 
negative in the Plains and Central States as far south as 
Arkansas. Maximum changes of 10° F. were reported 
near the Canadian border. Partially as a result of the 
new trough development, temperatures increased in the 
South, the Northwest, and the East. An increase of 6° F. 
occurred in Virginia in advance of the trough. 

The decreases of temperature in the Far West are more 
puzzling. An increase in 700-mb. height anomaly (fig. 7) 
does not usually imply a decrease in temperature. Many 
factors could contribute to a cooling, but perhaps most 
important in this case was the direction of anomalous 
flow, which changed from southwest the first half of the 
month to northwest the last 15 days, as the ridge retro- 
graded from the Rocky Mountains to the west coast. 
However, even though there was considerable intra- 
monthly cooling, mean temperatures were still above 
normal in the Far West and Southwest for the period 
January 16-30. 


4. WEATHER IN THE UNITED STATES 
RELATED TO CIRCULATION 


TEMPERATURE 


Several injections of Arctic air into the United Siates 
east of the Divide produced a cold regime in the !vast. 
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Ficure 6.—Mean 700-mb. contours (solid) and height departures 
from normal (dotted) (both in tens of feet) for January 16-30, 
1959. Note new trough which formed im central United States 
as the ridge in the West retrograded and the trough in the 
western Atlantic progressed. 
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Fieure 7.—Change of mean 700-mb. height departure from normal 


(tens of feet) from January 1-15 to January 16-30, 1959. This 
pattern accompanied considerable amplification of the long waves. 





This, together with warming in the West, resulted in a 
distribution of temperature which was meridional, a 
situation generally associated with great temperature 
extremes. In this month, for example, the warmest 
January of record was observed at Los Angeles and San 
Francisco, while Muskegon, Mich. and International 
Falls, Minn. were having the coldest (and snowiest at 
Muskegon) January of all time.! 

Extreme cold early in January accompanied a front 
which dropped daily temperatures as much as 50° and was 
responsible for a weekly average temperature of 20° 


—_—_—— 


' Fron al Climatological Data report of each city, January 1959. 





below normal at Amarillo, Tex. [8]. 
to the east coast, warming occurred throughout the West. 
An example of the magnitude of that change was recorded 
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Figure 8.—Temperature departure from normal (° F.) for (A) 
January 1-15, 1959, (B) January 16-31, 1959. Shaded areas are 
normal or above. In (A) warming east of Divide in Plains States 
was related to foehn action. In (B) note cooling in North Central 
and Northern Plains States. (C) Change in temperature depar- 
ture from normal (° F.) from January 1-15 to January 16-31, 
1959. Shaded areas represent positive changes (warming). 
(From [8].) 


As cold air spread 
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FicurE 9.—Percentage of normal precipitation for January 1959. 
Heaviest amounts were found near the mean 700-mb. wind 


(From [8].) 


maximum and cyclone track. 


at Great Falls, Mont. where a minimum temperature of 
—29° on January 3 was followed by a maximum of 62° 
on the 12th." By midmonth, warming progressed to the 
Plains States, while the Eastern States had become quite 
cold. During the week ending January 25, cold air again 
plunged to the Southern Plains States, following a fast- 
moving storm which became well organized in New 
Mexico on January 20. Average temperatures in Okla- 
homa, which were 6° above normal the preceding week, 
dropped to 6° below normal, and Birmingham, Ala. 
recorded a 52° decrease in temperature in 18 hours on the 
2ist. The reverse change in temperature was reported 
in the East from the Carolinas northward and resulted 
in a spectacular January thaw in the manner described 
by Wahl [9]. 

January was not especially noteworthy regarding 
temperature records, but there were additional features 
of interest. A comparison of figures 3 and 4 shows that 
the isopleth of 0° temperature departure from normal is 
almost identical with the zero isopleth of the average 
thickness DN. One also sees that below normal thick- 
nesses almost everywhere match below normal temper- 
atures. It should also be noted that the height anomalies 
of figure 1 have a similar relationship to the temperature, 
although less well defined. 

The 30-day mean of sea level pressure (fig. 2) also 
helps explain the temperature distribution of figure 3. 
The ridge from the Beaufort Sea to Florida and the flow 
from the District of Mackenzie to the Atlantic Coast 
States defines the source of cold air and its transport 
southeastward. The warm High over the Great Basin 
and the onshore flow into the West are indicative of the 
above normal temperatures observed there. It is of some 
historical significance that California coastal areas have 
escaped the below normal category since early April 
1958, except for two consecutive 5-day periods in mid- 
November. 
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Figure 10.— Mean sea level pressure departure from normal (milli- 
bars) and principal cyclone track (heavy dashed arrows). There 
was a close correspondence between the position of the storm 
track, the mean 700-mb. wind maximum, and area of heavy 
precipitation. 


PRECIPITATION 


Total precipitation accumulation in the United States, 
in terms of percentage of monthly normal, are shown in 
figure 9. The most prominent feature is the band of 
above normal precipitation which extended from the 
Pacific Northwest (where Spokane had its wettest Janu- 
ary), across Kansas, through the Ohio Valley, and into 
New England. Disastrous flooding occurred in Ohio, 
western New York (where Buffalo had its wettest Janu- 
ary), western Pennsylvania, and along the Wabash River 
system in Indiana. Rainfall of record proportions and 
melting snow associated with the storm of January 20-23 
were responsible for the flood loss of many lives and 
property loss of millions of dollars. 

It is generally harder to relate precipitation to the 
observed circulation or its anomalous flow than temper- 
ature. In this January, however, the task is relatively 
simple. In figure 10, the 30-day mean sea level pressure 
departures, one can see a corridor of negative values 
which closely resembles the configuration of the area 
of heavy precipitation. This was related to the principal 
storm track during January, represented by the heavy 
dashed arrow. The storm track is based on an analysis 
of frequencies of cyclones in a grid of equal area boxes 
and on inspection of tracks of daily centers of cyclones 
at sea level (Chart X of [5]). Note also the proximity 
of the 700-mb. mean wind speed maximum (in fig. }) 
to the cyclone track and axis of heavy precipitation. 
Compared with normal tracks and mean frequencies of 
cyclones [10], those in January were displaced southward 
5°-10° of latitude, presumably as a consequence of high- 
latitude blocking and the resultant depression of the 
westerlies. 
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Referring again to figure 9, other areas of above normal 
precipiiation are shown in California, where Mt. Shasta 
had the wettest January since 1916. There was a second- 
ary streak from northern Texas to South Carolina. The 
large area Which had less than 50 percent of normal 
precipitation included Dallas, Tex., which reported the 
driest January of record, and Albuquerque, N. Mex. 
which reported its sunniest. Dry weather in this area 
was related to northwesterly flow aloft south of the 
principal storm track. 
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Weather Note 


TORNADOES OF JANUARY 21, 1959--A FEATURE OF A WEATHER SINGULARITY? 


ROBERT R. DICKSON 
U.S. Weather Bureau, Nashville, Tenn. 


(Manuscript received February 2, 1959] 


ABSTRACT 


The tornadoes of January 21, 1959, in Tennessee and neighboring States are considered in relation to the New 


England January thaw singularity described by Wahl [1]. 
It is shown that at the time of the singularity there coexist 


to occur on the average on January 20-22 at Nashville. 


This singularity, in the form of a warm spell, is shown 


on the average in the Tennessee area certain conditions favorable for the formation of severe storms. These include 
a tongue of warm, moist air at the surface, a wind shift from southerlies during the warm period to cool northwester- 
lies immediately afterward, a 500-mb. trough to the west with southwesterly winds and contour inflection point 
over the Tennessee area, and the presence of a jet stream aloft. Review of past records reveals that tornadoes in 
the Arkansas, Alabama, Mississippi, Tennessee, and Kentucky area have occurred more often during the time of 
the singularity, January 20-22, than at any other time of the month. 


1. INTRODUCTION 


On January 18, 1959, surface winds from the Gulf Coast north- 
ward to the Great Lakes shifted to southerly quadrants, heralding 
the end of the existing cold spell. From Kentucky southward this 
warming continued through the morning of January 21, culmi- 
nating in temperatures of 70° F. or higher over most of the South- 
east. The surface weather system responsible for the final north- 
ward surge of warm, moist air was an intense open wave on the 
polar front which moved rapidly from central Texas to near Detroit 
by 0000 cmt, January 22. As the associated cold front swept 
across the Southeast displacing maritime tropical air, a series of 
severe thunderstorms and tornadoes broke out. Preliminary 
reports placed tornado activity in Mississippi, Alabama, Ten- 
nessee, and Kentucky. As the tornado-spawning Low moved 
northeastwvrd, 60° F. temperatures reached northward to Boston 
prior to the cold front passage near noon on January 22. 


2. THE NEW ENGLAND JANUARY THAW SINGULARITY 


The advent of 60° F. temperatures in the Boston area on Jan- 
uary 22, 1959, occurred at precisely the period (January 20-23) of 
the New England January thaw singularity described by Wahl [1]. 
Compared with the normal January daily maximum temperature 
at Boston of 36.5° F., this January 22 warming constituted an 
exceedingly well-developed example of the singularity. 

Wahl discussed the validity of the January thaw singularity for 
the eastern portion of the United States, using data for the Boston, 
New York, Washington, D.C., and Columbia, Mo., areas. Evi- 
dence of this singularity in Tennessee may be found in figure la. 
Mean daily temperature at Nashville, averaged over the period 
1871-1950 [2], shows an increase to 41.9° F. on January 21—nearly 
3° warmer than the average temperature 2 days before or 3 days 
after that date. Thus the Tennessee counterpart of the New 
England January thaw singularity is established as occurring on 
January 20-22. 

3. SYNOPTIC IMPLICATIONS 


Wahl [1] illustrates the remarkable extent to which the singu- 
larity shows up in 40-year average sea level pressure patterns. His 
figure 4b shows the existence of warming southerly winds on Jan- 


uary 20 from the Gulf of Mexico to New England, accompanying 
the warm-spell singularity on January 20-22 at Nashville and 
January 20-23 in New England. By January 27, as indicated in 
his figure 4c, the warm air over the eastern half of the country has 
been swept away by a cold High from the northwest. Figures 1b 
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Figure 1.—Data for Nashville, Tenn. (a) Daily average tempera- 
ture, 1871-1950; (b) Prevailing wind direction, 1931-1950; 
(c) Frequency of days with measurable precipitation, 1911-1950. 
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and 1c serve to document this sequence of events further. Pre- 
yailing vind at Nashville [3] is from the south on January 20-22 
giving «ay to northwesterly winds on January 23, thus suggesting 
a cold ‘ront passage. It is of interest to note that on January 21 
prevailing winds at Nashville are from the quadrant SW-SE 70 


percell of the time. Figure le shows a secondary maximum in 
the frequeney of days with measurable precipitation at Nashville 
(3] on January 21 and 22 accompanying the wind shift described 
aboves 

To further assay the potential of the synoptic situation associated 
with the January 20-22 warm-spell singularity, the concurrent 
upper-level circulation was investigated. Figures 2 and 3 are 500- 
mb. mean charts and their departures from normal for January 20 
and January 27 respectively for 1948-1957 (1955 missing). The 
500-mb. trough on January 20 stretches from western Texas north- 
ward to the Great Lakes, a position near the January normal. How- 
ever, the anomaly lines show that the southern portion of the trough 
is deeper than normal and winds over the eastern part of the country 
are more southerly than the January normal. Attention is also 
directed to the location of the 18,400-ft. contour on January 20 in 
the Arkansas, Kentucky, Tennessee area. As pointed out by 
Fletcher [4] the axis of maximum wind at 500 mb. is most frequently 
associated with the 18,400-ft. contour and is within the contour 
range 18,200—18,800 ft. 91.8 percent of the time. Since experience 
verifies a reasonably close identity between the location of the axis 
of maximum wind at 500 mb. and the jet stream core over the 
eastern United States, the jet stream on January 20 (fig. 2) can be 
considered to cross the Arkansas-Kentucky-Tennessee area near the 
18,400-foot contour. 

Figure 3 illustrates the eastward progression, shear, and decrease 
in intensity of the mean trough in central United States by January 
27. Height change from the January 20 average map of figure 2 
to the January 27 average map of figure 3 is shown in figure 4. 
The difference between these average daily 500-mb. charts one 
week apart is striking. 

As shown by the long-term average circulation at both sea level 
and 500 mb., some of the necessary conditions for the generation 
of severe storms in the general area exist at the time of the January 
20-22 warm-spell singularity in Tennessee. A surface tongue of 
warm (and by inference moist) air, introduced by persistent south 
winds, stretches from the Gulf through Tennessee to New England. 
This tongue of warm, moist air is subsequently displaced by cold 
air which brings a notable temperature drop and shifts surface 
winds from the south to the northwest in the Tennessee area. Aloft, 
a 500-mb. trough, deeper than normal in southern portions, lies 
immediately to the west, with southwesterly flow and contour in- 
flection point over the Tennessee area. Furthermore, the axis of 
the jet stream, as deduced from 500-mb. contours, traverses the area. 


4. TORNADO DISTRIBUTION IN JANUARY 


A search of past records from 1931 through 1959 for the Alabama, 
Mississippi, Arkansas, Tennessee, Kentucky area [5] reveals a total 
of 28 tornado days during the month of January. Of these, six 
(21 percent) occurred from January 20-22 in association with the 
singularity under consideration. A similar concentration of tornado 
occurrence is found in only one other January period, January 29-31, 
when six tornado days have been observed. The number of tor- 
nadoes on any January day for the 1931-1958 period reaches a 
maximum of 12 on January 22 and does not exceed 5 on any other 
January day. Furthermore, preliminary, unverified reports indi- 
cate that 11 tornadoes occurred on January 21, 1959, to give a 
preponderance of January tornadoes during the January 20-22 
warm-spell singularity. 


5. CONCLUSIONS 


Long-term average temperature data for Nashville give evidence 
lecnnessee area of a warm-spell singularity on January 20-22— 
the Tennessee counterpart of the New England January 
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FicurE 2.—500-mb. mean map for January 20, 1948-1957 (less 
1955). Height contours (100’s of feet) are solid lines, departures 
from normal (10’s of feet) are dashed lines. 
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Figure 3.—500-mb. mean map for January 27, 1948-1957 (less 
1955). Height contours (100’s of feet) are solid lines, departures 


from normal (10's of feet) are dashed lines. 











Figure 4.—500-mb. height change (10's of feet) from January 20 
mean map (fig. 2) to that of January 27 (fig. 3). 


thaw singularity described by Wahl [1]. Examination of average 
circulation features at sea level and the 500-mb. level reveals the 
existence of certain conditions favorable for the development of 
severe storms in the Tennessee area at the time of the January 
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20-22 warm-spell singularity. These include the presence of a 
warm, moist tongue of air in the lower atmosphere, passage of a 
cold front as suggested by a surface wind shift from persistent 
southerly winds prevalent during the warm spell to cooling north- 
westerly winds, a 500-mb. trough to the west with southwesterly 
winds and contour inflection point over the Tennessee area, and 
finally the presence of a jet stream aloft. It is suggested that in 
years when the singularity is especially well developed, as in Janu- 
ary 1959, the singularity and its attendant circulation patterns at 
sea level and aloft provide the ingredients for the generation of 
severe storms in Tennessee and vicinity. 

The results of this study are of value primarily as another bit of 
evidence dealing with the validity of a January 20—22 singularity in 
weather elements. Although the effect shows up more or less con- 
sistently in the long-term averages, it is not reliable enough to de- 
pend on every year. Because other sources of large variability are 
present, the singularity has limited value as a forecasting tool. It 
serves primarily to focus attention on the potential of the period. 

This study was limited to an examination of the January 20-22 
singularity in data for the Tennessee area in relation to the regional 
circulation singularity described by Wahl [1]. How the results fit 
into the worldwide pattern of January singularities suggested by 
several recent studies is a tempting question for further speculation. 
It is interesting to note, for example, that Brier’s [6] data on hemi- 
spheric fluctuations in the meridional exchange of air at 50° N. 
latitude, Bowen’s [7] rainfall data for a large number of stations in 
the Northern and Southern Hemispheres, Bigg’s [8] high-level cloud 
data at Australian stations, and Kline and Brier’s [9] freezing nuclei 
counts at Washington, D.C., also were peaked around January 


20-24. 
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NOTICE 


Effective with this issue (January 1959) the 17 climatological charts 
ordinarily inserted at the end of each issue will be discontinued in the Monthly 


Weather Review. 


They will continue to be printed (in black and white) in the 


Climatological Data, National Summary which is also issued monthly. 

Paid subscribers to the Review, who do not already receive the Climato- 
logical Data, National Summary, will be added to the CDNS mailing list for 
six months or until their present Review subscriptions expire, whichever is earlier. 
Cooperators who need the charts will be placed on the mailing list for the 
Climatological Data, National Summary upon request to Chief, U.S. Weather 


Bureau, Washington 25, D.C. 
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